NON-ANALYTICITY AND THE VAN DER WAALS LIMIT « 

S. Friedli and C.-E. Pfister 

Abstract: We study the analyticity properties of the free energy /^{m) of the 
Kac model at points of first order phase transition, in the van der Waals hmit 
7 \ 0. We show that there exists an inverse temperature /?o and 70 > such that 
for all (3 > I3q and for all 7 G (0,70), f-yim) has no analytic continuation along 
the path m \ m* (m* denotes spontaneous magnetization). The proof consists 
in studying high order derivatives of the pressure p-yih), which is related to the 
free energy f^{m) by a Legendre transform. 
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1. Introduction 

The first equation of state giving precise predictions on tlie liquid-vapor equi- 
librium at low temperature was given by van der Waals |vdWj : 

(^+^)("-^)=^^- (^-^^ 
This equation follows from the hypothesis that the molecules interact via 1) a 
short range hard core repulsion, due to the assumption that molecules are ex- 
tended in space, 2) an attractive potential, whose range is assumed to be com- 
parable to the size of the system. Nowadays, such an approximation is called 
a mean field approximation. As well known, there exists a critical temperature 
Tf. = Tc{a, b) such that for T < T^., -^p > for some values of v, which implies 
thermodynamic instability. On physical and geometrical grounds, the graph of 
the pressure was modified by Maxwell who replaced p{v), on a suitably chosen 
interval [1;^,^;^], by a flat horizontal segment (the "equal area rule"). The new 



MCp{v) 
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Figure 1: The equation of state modified by Maxwell and the analytic continuation at 
the condensation point. 

function obtained, written MCp(t'), describes precisely what is observed in the 
laboratory: vi is called the evaporation point and Vg is the condensation point 
(see Figure H}. 

A particularity of this scenario is that MCp can be continued analytically along 
the paths v y vi and f \ w^: the liquid and gas branches can be joined analyti- 
cally by a single function, which is nothing but the original isotherm p given in 
The pressure obtained by analytic continuation was originally considered 
as the pressure of a meta-stable state (see Figure [T]). For instance, the meta- 
stable state obtained by analytic continuation along the path v \ Vg is called a 
super-saturated vapor. 

Much later, Kac, Uhlenbeck and Hemmer |KUHj showed how the Maxwell 
construction could be rigorously justified for a one dimensional model, from first 
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principles of statistical mechanics, using a double limiting process: if the range of 
interaction diverges after the thermodynamic limit, then convexity is preserved 
and the free energy converges to the convex envelope of mean field theory. Later 
this was generalized and extended to higher dimensions by Lebowitz and Penrose 
|LPj . From the point of view of analyticity, these results imply, as in the theory 
of van der Waals, that the free energy can be continued analytically across con- 
densation/evaporation points. 

In the mean time, arguments were given, saying that when the range of interac- 
tion is finite, the free energy might have some singularities that forbid analytic 
continuation across the transition points. In jFj and |Llj . Fisher and Langer 
analyzed in details simple models to illustrate this phenomenon, but it was not 
until the seminal work of Isakov jTT] that this was shown for the Ising model. 

An important issue is thus to understand how the breakdown of analyticity 
at a first order phase transition point relates to the range of interaction. Since 
Kac potentials give a way of interpolating finite range systems and mean field, it 
seems an interesting problem to study the dependence on the scaling parameter 
7 of the analyticity properties of the Kac model at low temperature. The aim 
of this work is to show that as long as this parameter is strictly positive, i.e. as 
long as the range of interaction is finite, the free energy has no analytic contin- 
uation at first order phase transition points. Only after the van der Waals limit 
(7 \ 0) does the free energy have analytic continuation. This result answers 
a question raised by Joel Lebowitz at a conference devoted to Kac potentials. 
In-homogeneous Random Systems, held in Paris, January 2001. 

In Section 11.11 we remind the main properties of the free energy for mean field 
and Kac potentials in the case of Ising spins. In Section 11.21 we state our main 
results and give the strategy of the proof. 

1.1. Mean Field and Kac Potentials. We consider the lattice Z"^, d>2, with 
a distance d{x, y) = \\x — y\\, where 

||a:|| := max . (1-2) 

1=1,. ..,d 

This distance will also be used for points of M.^. The letter A will always denote a 
finite subset of Z'^. At each site z G Z"^ lives a spin cTj G {±1}- The configuration 
space is r2 = {±1}'^ . For any set A, Q\ = {±1}^. Our notations are often 
inspired by those of Presutti |Prj . 

Mean Field. In a mean field model, the interactions ignore the spatial positions 
of the spins, and the hamiltonian in a volume A containing sites is {a G Qa) 

0):=-l Yl (1-3) 

{«J}CA 

As is well known, the free energy can be easily computed. For m G [—1, +1], 

/a/fM = -^m^ - i/(m) , (1.4) 
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where 

^, . 1— m 1 — m 1 + m 1+m 

/(m) := ^loS^ 2~ ^^^~' ^ 

When /5 < 1 /a/f is strictly convex, but when /5 > 1, /mf has two minima 
at ±m*(/3), where m*(/?) is the positive non-trivial solution of m = tanh(/3m). 
/3c := 1 is the critical temperature of mean field theory. As in van der Waals 
theory, f^p is non convex when (] > (3^-, in contradiction with thermodynamic 
stability. 

Kac Potentials. Kac potentials are defined as follows. Consider J -.W^ ^ M"*" 
supported by the unit cube {y G W'- : \\y\\ < 1} = [—1, +1]'^ such that the overall 
strength equals unity, i.e. 



J{x)dx = 1 . (1.6) 

Let 7 G (0, 1) be the scaling parameter. Define -.U^ IR+ as follows: 

J^[x) := c^-f'^J{'yx) , (1.7) 

where is defined so that 

^J,(x) = l. (1.8) 

It is easy to see that p.6p implies lim^x^oc-y = 1. Since J^{x) = if > 
we call R := 7"^ the range of the interaction. Unless stated explicitly, R will 
always denote the range of interaction, i.e. 7"^. For a finite A, a G Qa, the Kac 
hamiltonian is defined by 

Hl{^) = - Yl M^-J)(^^^J-hJ2^^^ (1-9) 
where /i G M is the magnetic field. The magnetization in A is 



'ieA 



and takes values in a set xa C [—1, +1]. The canonical partition function is defined 
by (/3 > is the inverse temperature, m G xa): 

Z{A,m)= Yl exp ( - /9iJ°(a)) . (1.11) 

'nA(o-A)=m 

The free energy density is, for m G [—1, +1], 

f,{m) = - hm -i-logZ(A,m(A)), (1.12) 

where the thermodynamic limit A /" Z*^ is along a sequence of cubes, and the 
sequence m(A) is such that m(A) — > m. The function exists and is convex. 
The Theorem of Lebowitz- Penrose |LPj gives a closed form for the free energy in 
the van der Waals limit 7 \ 0. For a function f{x), let CE f{x) denote its convex 
envelope. 



5 



Theorem 1.1. |LPj For any (3 > 0, m e [-1, +1], 

/o(m) := lim/^(m) = CE fMF{m) . (1.13) 

When P > 1, the graph of fo{m) is thus horizontal between —m*{(3) and +m*{j3), 
giving a rigorous justification of the MaxweU construction (see Figure E]). 



-1 



Figure 2: The free energy fo{m) when (3 > 1. The dotted Hne is the analytic continu- 
ation provided by fMpini). 



From the point of view of analyticity, we have 

Corollary 1.1. When (3 > 1, fo is analytic everywhere except at ±m*(/3), and 
has analytic continuations along the (real) paths m —m*{(3), m \ +m*{f3). 
The unique analytic continuation is given by the mean field free energy /mf- 

That is: after the van der Waals limit, all the analyticity properties of the free 
energy are known explicitely. There exists no formula for f^ when 7 > 0, and it 
was not shown, until the papers of Cassandro and Presutti |CPj and Bovier and 
Zahradm'k |BZlj . that the system exhibits a first order phase transition before 



reaching the mean field regime: for all > 1, the graph of f^yim) already has a 
plateau [—m*{j3, 7), +m*{/3, 7)] when 7 is small enough. In this sense, one can say 
that mean field, together with the Maxwell construction, is a good approximation 
to long but finite range interactions (and vice versa). Our purpose is to show that 
from the point of view of analyticity, the situation is very different. 

1.2. Obstruction for 7 > 0; Main Results. Our results hold for Kac poten- 
tials for which Lemmas 12.11 and 12.81 hold, but we believe them to be true for any 
ferromagnetic potential satisfying ()1.6p . For the sake of simplicity, we focus on a 
particular potential, i.e. on the step function 

J(a:) :=2-'^l|||,ll<i|(x), (1.14) 

In this setting, our main result for the free energy density is the following: 

Theorem 1.2. There exists (3q and 70 > such that for all (3 > (3q, 7 G (0, 70), 
is analytic everywhere except at ±m*{f3,'j), but has no analytic continuation 
along the paths m y —m*{P,j), m \ -|-m*(/?,7). 

This result is in favor of the original ideas of Fisher and Langer, saying that 
finiteness of the range of interaction is responsible for absence of analytic con- 
tinuation. In particular it excludes the possibility of obtaining the free energy 
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by a Maxwell construction: when 7 > the phases + and — cannot be joined 
analytically. 

The proof of Theorem 11.21 will be done by working in the more appropriate 
grand canonical ensemble (in the lattice gas terminology), in which the constraint 
on the magnetization is replaced by a magnetic field. Let 

Z(A)= exp ( - /3i/;j(aA)) . (1.15) 

Define the pressure density by 

p^{h) := lim p^^A^h) , where p-y,A{h) = --^logZ(A) . (1.16) 

The free energy and pressure densities are related by a Legendre transform: 

/y(m) = sup(/im — p-y{h)) . (1-17) 

See for instance |Prj for a proof of this property. The analytic properties of 
at ±m*(/3,7) will be obtained from those of p^ at h = 0. By the Theorem of 
Lee and Yang |LYj . p^ is analytic outside the imaginary axis. The main result 
of the paper is the following characterization of the analyticity properties of the 
pressure at h = 0. 

Theorem 1.3. There exists (3q, 70 > and a constant > such that for all 
f3 > Po, 1 ^ (0,7o); the following holds: 

1) The directional derivatives p^^''~{0) exists for all A; G N, i.e. p^ is C°° at 
h = 0. Moreover, there exists a constant C+ > such that for all /c e N, 

sup '^(/i)! < (C+7^/3-^)'/c!^ + C^kl . (1.18) 

0<Rc/i<e 

2) The pressure has no analytic continuation at h = 0. More precisely, there 
exists > and an unbounded increasing sequence of integers ki,k2, ■ ■ ■ such 
that for all k G {A;i, k2, ■ ■ ■ }, 

'^(0)1 > (C_7^/5-^)''A;!^ - C^k\ . (1.19) 

The lower bound ()1.19|1 becomes irrelevant when 7 \ 0. Moreover, we should 
mention that each integer ki depends on 7 and P, with lim^\o h = +00: informa- 
tion about non- analyticity is lost in the van der Waals limit. Since we know from 
the Lebowitz-Penrose Theorem that p^ converges, when 7 \ 0, to a function 
that is is analytic at /i = 0, it is worthwhile considering the low order derivatives 
of p^/. Considering the upper bound ()1.18p . it easy to show the 

Corollary 1.2. There exists C = C{(3) such that for small values of k, i.e. for 

k < 7"°', we have the upper hound 

sup \p^J;^'*-{h)\<C^k\, (1.20) 

0<Re/i<e 
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Figure 3: The derivatives of the pressure at h = 0, when 7 > 0. The first ones (fc < 
7"'') behave hke those of an analytic function, but non-analyticity always dominates 
for large k. 



This shows that a close inspection of the derivatives of the pressure allows to 
detect how analyticity starts to manifest when 7 approaches 0. These different 
behaviours are illustrated on Figure El 

To show Theorem II. we first construct the phase diagram of the Kac model 
with a complex magnetic field, at low temperatures, 7 small. Then, we adapt the 
technique of Isakov to obtain lower bounds on the derivatives of the pressure in 
a finite volume. These two essential steps deserve a few comments. 

(1) Phase diagrams of lattice systems can be studied in the general frame- 
work of Pirogov-Sinai Theory ( |PSj . |Zlp . which applies when the system 
under consideration has a finite number of ground states, and for which 
the unperturbed hamiltonian satisfies the Peierls condition. In our case, 
the Kac potential has two ground states which are the pure + and pure — 
configurations, but the Peierls constant (computed with respect to these 
two ground states) goes to zero when 7 \ since in the van der Waals 
limit, the interaction between two arbitrary spins vanishes. Therefore, a 
direct application of Pirogov-Sinai Theory would lead to a range of tem- 
perature shrinking to zero in the van der Waals limit. 
Recently, Bovier and Zahradnik |BZ2j proposed a systematic method to 
study spin systems with long but finite range interactions. Their technique 
allows to study, for instance, the Kac model with a complex magnetic field, 
in a range of temperature that is uniform in 7. In their approach, the 
ground states of Pirogov-Sinai Theory are replaced by restricted phases, 
i.e. by sets of configurations. In the +-restricted phase, for example, all 
the points are -I— correct, i.e. their 7" ^-neighbourhood contains a majority 
of spins -|-. When a point is in neither of the restricted phases, it is in 
the support of a contour F, and it can then be shown that the contours 
defined in this way satisfy the Peierls condition with a Peierls constant 
p that is uniform in 7: ||F|| > p|F| where ||F|| is the surface energy of 
F. In Section El we show that a polymer representation can be obtained 
for the restricted phases, and that their corresponding free energies be- 
have analytically aX h = 0. The full phase diagram is then completed 
in Section jH we give precise domains in which the partition function 
can be exponentiated. These domains are made optimal by introducing 
special isoperimetric constants associated to contours (see the discussion 
hereafter, and ()2.47|) ). Complications arise from the fact that polymers 
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of the restricted phases induce interactions among contours. Besides the 
definition of the restricted ensembles, our analysis of the phase diagram is 
independent of the paper |BZ2j . In a different setting, restricted ensembles 



were also studied in |BSj . |DSj and |BKLj . 

(2) To implement the mechanism used by Isakov, we consider the pressure 
in a finite box A, with a pure +-boundary condition. By introducing 
an order among the contours inside A, the pressure can be written as a 
finite sum: 

P\ \ PI I rec+(A) 

where ^^(A) is the restricted partition function and C''~(A) is the family 
of all contours of type + in A. With the analysis of Sections El and El 
the derivatives of the functions u~}^{r) can be estimated using a stationary 
phase analysis. When A is sufficiently large, the contributions to p^^\{0) 
are the following: because it is analytic, the restricted phase contributes a 
factor C^k\. Then, a class of contours called /c-large gives a contribution 

d 

of order kl ^--^ . The rest of the contours is shown to have a negligible 
contribution in comparison of the A;-large ones. This gives a lower bound 

\p^!a\0)\ > (C„7^r^)'A;!^ - C^k\ . (1.22) 

In the last step of the proof we show that limA p^'l'^'* (0) = P7^^'^(0), and 
so p.22|) extends to the thermodynamic limit A Z"^, which gives ()1.19|) . 

Before going further, we make an important remark. In jl^, Isakov proved 
Theorem 11.31 for the Ising model. An attempt was then made, in a second paper 
|I2] . to extend the method to any two phase model for which the Peierls condi- 
tion holds. Unfortunately, this extension could only be done under two additional 
assumptions which we briefiy describe. Associate to each phase a discrete isoperi- 
metric problem of the following type: let V{T) denote the volume of the contour 
r (of a given type) and ||r|| its surface energy. For G N, consider the problem: 

'Find the best constant C{N) such that ^ < C{N)V{T)-d 
for all contour F with V(V) < N . 

The assumptions of Isakov are then that in the limit oo, 1) the asymptotic 

behaviour of the constant C{N) is the same for the two phases, 2) there exist 
maximizers of arbitrary large size. 

Clearly, these assumptions are satisfied by the Ising model, for which ||F|| = |F| 
(the number of dual bonds on the dual lattice) and the maximizers are always 
given by cubes, i.e. C{N) = (2d)~^ for all A^. But for a model with no symmetry 
or with interactions that are more complicated than nearest neighbours, these 
assumptions can be very hard to check. The problem comes from the fact that 
the surface energy ||F|| depends on the detailed structure of the hamiltonian. In 
our case, symmetry reduces the difficulty to the existence of large maximizers. 
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We will see that the construction of the phase diagram can be done when the 
isoperimetric problem is formulated as follows: 

Find the best constant K{N) such that ^ < K{N)V{T)Tt 
for all contour F with V^(F) > . 

By formulating the problem in this way, the existence of large maximizers is 
immediate, and we avoid the necessity of solving the isoperimetric problem ex- 
plicitly. 

It was actually shown in |FPj that the two assumptions of Isakov can be swept 
out, and that the result of [l2] can be extended to the whole class of two phase 
models treated generally in Pirogov-Sinai theory, the only necessary ingredient 
for non-analyticity being the Peierls condition. The general theorem of |FPj ap- 
plies to the Kac model but with some restriction (3 > Po{l) where Polj) diverges 
when 7 \ 0. In the present paper we study the van der Waals limit at fixed jS. 

The description of the model in terms of contours and the verification of the 
Peierls condition for ||F|| will be done in Section|21 Section|21is entirely devoted 
to the study of restricted phases and to their analyticity properties, adapting the 
technique of |BZ2j . Section E] is the construction of the phase diagram in the 
complex plane of the magnetic field. Section contains the proofs of our main 
results. In Appendix^ we give the details of the stationary phase analysis for the 
study of the derivatives of the functions u~l{r), and Appendix iBl contains basic 
definitions for the cluster expansion technique. 

Conventions: we will often use the norm := sup^g^, When G is a 

graph we denote by V{G) its set of vertices and by E{G) its set of edges. 

Acknowledgments : We wish to thank Anton Bovier and Milos Zahradm'k for many 
useful discussions concerning |BZ2j . 



2. Contour Description 

For the description of configurations in terms of contours, we use the notion of 
correct /incorrect point introduced by Bovier and Zahradnik in |BZ2j . There are 



two major requirements for the way in which contours should be defined. 

(1) They are defined on a coarse-grained scale, and a Peierls condition must 
hold for the surface energy of each contour, with a Peierls constant that 
is uniform in 7. See Proposition 12.21 

(2) Outside contours, a partial re-summation over configurations will lead to 
restricted phases. To obtain convergent expansions for these phases, care 
must be taken in the definition of contours. See the parameter 6 in (j2.16p . 

2.1. Definition of Contours. We introduce some more notations. We have 
d{x,A) = M{d{x,y) : y e A}. For > 1, define the box Bn{x) := {y e Z'^ : 
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d{x,y) < N}, and B^{x) := Bi^{x)\{x}. The A^-neighbourhood of A is 

[A]^:= IJS^(x), (2.1) 

and the boundaries 

d+A = {x e A' : d{x,A) < N}, (2.2) 
d'A = {xeA: d{x, A") < N} . (2.3) 

A set A is A^-connected if for all x, y G A there exists a sequence Xi,X2, ■ ■ ■ , Xn-i, Xn 
with xi = X, Xn = y, Xi G A, and (i(xj, Xj+i) < A^. If a\ G I^a, Vh" ^ ^A'^; we 
define the concatenation a^rj^c g f2 in the usual way: 



(-a^ao.= r^;^ (2.4) 




We often use the symbol # to denote either of the symbols + or — , or the constant 
configuration taking the value # at each site of Z*^. We define 

0ii(o-i, aj) := -^J^{i - j){(ri(rj - 1) , (2.5) 

Let (pij := (pij{+,—). The overall interaction strength is the upper bound on the 
energy of interaction of a single spin with the rest of the system, and equals 

5^0., = 5^ J,(z-j) = l. (2.6) 

Relevant functions for the study of nearly constant spin regions are the following 
(they will appear naturally later when reformulating the hamiltonian) : 

wfji'^i^ ^j) := <Piji(^i, CTj) - o-j) - (pijiai, #) . (2.7) 

Notice that w*{i^,aj) = w*{ai,if) = 0. Let S G (0, 1), a E Q. With regard to 
the step function J defined in (|1.14|1 . we define a point i to be {6, +)-correct for 
a if 

|i?^(^)n{J:a, = -l}|<||SK(^)|. (2.8) 

That is, the /^-neighbourhood of a (<5, +)-correct point contains a majority of 
+ spins. Although we will always consider the step function, it is often easier 
to formulate proofs with the help of the functions wf^, since they will appear 
naturally later in the re-formulation of the hamiltonian. We thus define the 
notion of correct /incorrect point in the general case. 

Definition 2.1. Let 6 G (0, 1), a e n, i e TJ^ . 

(1) % IS (5, +)-correct for o ifT^j-.j^i l^tji-, < ^■ 

(2) i IS (5, -)-correct for a ifY^j-.j^i < ^■ 

(3) i is 5-correct for a if it is either (<5, +)- or (<5, —)-correct for a. 

(4) i is 5-incorrect for a if it is not 6-correct. 
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It is easy to see that this definition coincides with ()2.8p when J is the step 
function. 

The notion of correctness for a point i depends on the spins in the i?-neighbour- 
hood of i but neither on the value of a,, nor on the magnetic field. Notice that if 
5 = this notion of correct point essentially coincides with the one of Zahradnik 
in |Zlj . We first show that when 5 is small, regions of (5, +)- and (5, — )-correct 
points are distant. In particular, a point i cannot be at the same time (5, +)- and 
((5, — )-correct. 

Lemma 2.1. Let 5 e (0,2"'^), a ^Vt. Then 

1) Ifi is {6, +)- correct, the box Bji{i) contains either {6, +)- correct, or 6 -incorrect 
points (but no {5, —)-correct points). 

2) Ifi is (5, —)-correct, the box BR{i) contains either (5, —)-correct, or 5 -incorrect 
points (but no (5, +)-correct points). 

Proof. Suppose i is (5, +)-correct for a. Consider j G BR{i) and compute 

E l^7fc(+'^'^)l= E 20.^^ E 20,', (2.9) 

(Tfc= + 1 CTfc= + l 

Using the properties of the function J ^, we can exchange j and i and write 

^ 20,fc = ^ 20,fc = ^ 2(j),k - 20ifc (2.10) 

fcGB^(i)nB*(i) keB'^ij)nB'^{i) k^i k^B'^{j)nB'j,{i) 

fTfe=+i fffc=+i '7fc=+i fTfc=+i 

Using ()2.6|) and |i?_R(i) H BR{i)\ > 2^'^\BR{i)\, this last sum can be bounded by 

E (2-11) 

Then, since i is (<5, +)-correct for cr, 

E 20,fc = 2- E 20,, = 2- E |t/;+(-,a,)| >2-5. (2.12) 

0-fc= + l o-fe=-l 

We thus have the lower bound 



5^ \w-,{+,a,)\>2-5-2^^>5, (2.13) 

i.e. J cannot be {6, — )-correct for a, which finishes the proof. □ 

In the sequel we will always assume that 6 G (0,2""') is fixed. The cleaned 
configuration a G ^2 is defined as follows: 

-1 if i is (S, +)-correct for a , 
ai := { —1 if ^ is {5, — )-correct for a , (2-14) 
(jj if i is (5-incorrect for a . 



^At this point we use the particularity of the step function: (l)jk is constant on the intersection 
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For any set M C Z'^, we can always consider the partial cleaning a]\jWMc which 
coincides with cr on M and with a on M'^. In the sequel, the cleaning and partial 
cleaning are always done according to the original configuration a, with a fixed 
6. Notice that if a point i is, say, {6, +)-correct for cr, then the cleaning of a has 
the only effect, in the box Bf}{i), of changing — spins into + spins (and not + 
spins into — spins). This is a consequence of Lemma l2.ll We denote by /^(o") 
the set of (5-incorrect points of the configuration a. The important property of 
the cleaning operation is stated in the following lemma. 

Lemma 2.2. Let Mi C M2, 5' G [0,6]. Then Is'io-M.aA.q) C h'i^^Ah'^A'q) ■ 

Proof. Let z be a (5', +)-correct point of crAf2^M|- Using the fact that ctmiO'm^ 
and (JM^^aM^ coincide on Mi and M|, we decompose 

k:k^i k:k^i k:kj^i 

There are at most three possibilities for a point k of the last sum. 1) If A; is 
((5, +)-correct for a then ak = +1 and so \w^ki~^^k)\ = 0. 2) If k is 5-incorrect 
for a then Wk = (Tk = {o'M2^M§)k- 3) If k is (5, — )-correct for a then it is also 
(5, — )-correct for (Tm^o'm^- By Lemma EUl i is not (5, +)-correct for (Tm2'^m;^- 
This is a contradiction with the fact that i is (5', +)-correct for o'M'i'^Aq, so there 
are no such k. 

We can then bound the whole sum by 6'. This shows that i is {6', +)-correct for 
^Mi^Mf 5 and finishes the proof. □ 

Contours are defined on a coarse-grained scale. Consider the partition of Z'^ 
into disjoint cubes C^''^ of side length / G N, / > 2R, whose centers lie on the sites 
of a square sub-lattice of Z"'. We denote by the unique box of this partition 
containing the site i G 7i^. C^'-* will denote the family of all subsets of Z"' that are 
unions of boxes C^^-*. For any set A C Z^, consider the thickening (compare with 

m) 

{A}i := U C« . (2.15) 

In the sequel we always consider / such that / = uR, with u > 2. 
We will need to decouple contours from the rest of the system. Since interac- 
tions are of arbitrary large finite range, we follow |BZ2j and introduce a second 



parameter 6 G (0, S) . This new parameter is crucial; its importance will be seen 
later, for instance in the proof of the analyticity of the restricted phases. For 
each (T G r2 with |/5(cr)| < 00, consider the following set: 

S{a) := {M G C« : M D [/^(a)]^, M D [/KcTA/^A/0]if } • (2-16) 

First we show that S{a) is not empty. Consider Mq := {[Isi.'^)]R}i- If = then 
-^si^) — ^s{<^) = and any subset of Z'^ is in S{a). So we assume Mq 7^ 0. This 
gives £{a) since Mq G C«, Mq D [/^(a)]^ D [/^(a)]^ and Mq D [Tsia)]R D 
[Ig{crMo^M§)]R by Lemma We then show that S{a) is stable by intersection. 
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Suppose A, B E S{a). Then clearly An B D [h{.o-)]R and using again Lemma 

A D [Is{aAaAc)]R D [Isi(^AnB^{AnBr)]R , (2.17) 

B D [IsiaBaBc)]R D [/^(crAnBO^CAnB)^]^ > (2-18) 

which implies An B E ^(^)- The following set is thus well defined, and is the 
candidate for describing the contours of the configuration a: 

r{a):= Pi M. (2.19) 

Mee{a) 

By construction, I*{cr) is the smallest element of S{a). A first important property 
of /*(o") is the following, which will be essential to obtain the Peierls bound on 
the surface energy of contours. 

Lemma 2.3. There exists, in the 2R-neighbourhood of each box C^^^ C I*{o'), a 
point j G which is 5 -incorrect for the configuration ai*aj*c. 

Proof Let C^'' C /*((t). First, suppose /^(a) n [C^%r 7^ 0. Then each^j G 
/^(cr) n [C^^^]2R is 5-incorrect for a, and hence 5-incorrect for ai*aj*c, since 6 < S 
and a and aj*ai*c coincide on Bji{j). 

Suppose there exists a box C^') such that ^ [/^(o-)]/^ n [C^^'>]r = 0. If I^{aj*aj*c) n 
[C(')]2ij = 0, i.e. [/^(a/.a/*.)]/? n [C(')]/j = 0, then we define /' := /*\C7(') and 
show that /' G S{cr), a contradiction with the definition of /*. First, /' D [Is{<7)]r- 
Using Lemma /* D [/5(o-/.a/.c)]^ d [/5(o-//CT//c)]r. Since [/^(a/*a/«)]R H 
[C^% = 0, this implies /' D [/5(a7.a7,c)]^, i.e. /' G £'(a). □ 

When studying restricted phases, we will need to re-sum over the set of con- 
figurations that have the same set of contours, that is to consider, for a fixed a 
(we assume I*{cr) 7^ 0), 

Aia) := {a' : a'j,^,) = r(a') = r(a)} . (2.20) 

It is important to have an explicit characterization of the set A{a). Let A*(cr) 
denote the set of points of I*{aY that are [5, #)-correct for a. By Lemma f2. II we 
have d{A~^ (a) , A~ (a)) > I, and we have the partition 

Z'^ = r(a)UA+(a)UA-(a). (2.21) 

We now show that the set A{(7) can be characterized explicitly by 

■D(cr) := {a : = 0-7.(0-), each i G [A*(cr)]ij is {S, #)-correct for a'^ . 

Proposition 2.1. ///*(a) ^ 0, then A{a) = V{a). 

Proof 1) Assume a' G A{a). Then /* = r{a) = r{a') D [/^(aOliJ. so that each 
i G [/*'^]_R is (5-correct for a'. Let A be a maximal connected component of [/*'^]_r. 
There exists i E A such that i E P, since we assumed P ^ ^. By Lemma l2.ll 
it suffices to show that i is (5, +)-correct for a if and only if it is (5, +)-correct 



^Hcre we use the fact that A n [B]2ii = if and only if [AJ^ n [5]^ = 0. 
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for cr'. Assume this is not the case, i.e. suppose i is {6, +)-correct for a and 
{6, — )-correct for a'. That is, 

^k+(-,(a,.a,..).)l= E Ki-,<^j)\<h (2.22) 
jy^i jeB*(i)n/* 

j^i jeB*(i)n/* 
Since z G /* we have ^ 

E k^-(+,(^/*^/-).)l< E l'"^.^-(+,+)l<2(l-2-'^). 

Therefore we get a contradiction, since, 

<25 + 2 E |%(+,(^T/.a/.Oi)l < 2<5 + 2(l -2"'^) < 2, (2-24) 
where we used the fact that 6 < 6 < 2^'^. 

2) Suppose a' G ^^(o"). Since a' coincides with a on I* (a) and all points of 
[I*{aY]R are (^-correct for a', we have /^(cr') = /^(cr). This gives I*{<j) D 
[^<5(c^)]i? = [h{cr')]R- Then, since (T/.(^)a7.((,)c = cr'/*(^)a'/.(<^)c, we have /*(cr) D 
[/5(cr/.(^)a/*(^)c)]ij = [/^(cr'/.^o^VHOli?- This implies J* (d) G £^(cr'), i.e. /*(cr') C 
I*{(y). Assume I*{a)\I*(a') ^ 0. Using the fact that a and cr' coincide on 
r{(r)\r{a'), we have cr/.(^/)CT/.(^/)c = o-/.(^/)^7.(^/)c. This gives, like before, 

r{a') D [/5(<7'/*(-')^'/*K)0]i?. = [l5i^iH-')^i'(-r)]R- With /*(a') D [/^(aOliJ = 
[/5(a)]ij, this implies J*(a') G ^(a), i.e. /*(a') D /*(a). So a' G ^(a). □ 

In particular, Proposition l2. ll implies that the cleaned configuration o"/*(cr)a/.(o-)c 
is an element of A{cr). 

Definition 2.2. . The connected components of I* (a) form the support of the 
contours of the configuration a, and are written supp Fi, . . . , supp r„. A contour 
is thus a couple T = (suppr,crr), where is the restriction of a to T. 
A family of contours {Fi, . . . ,r„} is admissible if there exists a configuration a 
such that {Fi, . . . , F^} are the contours of a ^. 

The fact that the contours are defined on a coarse-grained scale will be crucial 
when dealing with the entropy of contours, which we must control uniformly in 
7. Notice that two (distinct) contours are at distance at least / from each other. 
We will usually denote suppF also by F. Contours should always be considered 
together with their type and labels, which we are about do define. The following 
topological property is needed for the definition of labels. 



■^Hcrc we use a property of the step function, but this can be done for any Kac potential 
whose function J has the symmetry J{x) = J{y) when ||x|| = \\y\\. 

^Note that this configuration is not unique, unhke in the usual situation treated in Pirogov- 
Sinai Theory. 
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Lemma 2.4. Fix R > 1. Let B C be R-connected and bounded. Then d^A 
and d]^A are R-connected, where A is any maximal R-connected component of 
B" = Z\B. 

Proof. Let A be any maximal /^-connected component of B'^. Then A'^ is R- 
connected. Indeed, let x,y E A'', and consider a path xi = x, X2, ■ ■ ■ , Xn = y, 
Xj_|_i) < R. If Xi G A"^ for all i there is nothing to show. So suppose there 
exists I < i- < i+ < n such that {xi, . . . , C A'^, G A, G A, 

{xi^, Xi^+i, . . . ,Xn} C A'^. Since A is maximal, we have Xi_ G S, Xj^ G .B, and 
we can find a path from to Xi^ entirely contained in B, i.e. in A'^. 
We then show that dfA is -R-connected. Fix e > and consider the sets 

X = {x eR"^ : d{x,A) <^ + e} , (2.25) 

F = {yGM":%,A^)<| + e}. (2.26) 

Then X, Y are closed arc-wise connected subsets of M"', and X UY = R"'. By a 
Theorem of Kuratowski, X fl y is arc-wise connected ^. Let e' > and consider 
X, y G together with x,y E X HY such that (i(x, x) < |, (i(y, ?/) < |. Then 
consider any sequence Xi = x, . . . , x„ = y, Xj G X fl F, (i(xj, Xj+i) < e'. For each z 
we have d{xi, A) < ^ -\- e, d{xi, A^) < -| + e. This implies that each box BR_^_^(xi) 

contains at least one element x^ G diA, i.e. d{xi,x[) < -f + e. We have 

d{x[, x-_^J < rf(x-, Xi) + d{xi, Xi+i) + d{xi+i,x[^^) < R + 2e + e' . (2.27) 

If 2e -|- e' < 1, this shows that dfA is i?-connected, which implies that d~I^A is 
i?-connected. The same proof holds when d'^A is replaced by d]^A. □ 

Let r be a contour of cr, A a maximal /^-connected component of (supp F)'^. Let 
i G d]^A. By definition, i is [S, #)-correct for a for some # G {±1}- By Lemmas 
12.41 and 12.11 each i' G d]^A is (5, 7^)-correct for cr for the same value We call 
^ the label of the component A. There exists a unique unbounded component of 
F'^. The label of this component is called the type of the contour F. Let F be of 
type + (resp. — ). The union of all components of F'^ with label — (resp. -|-) is 
called the interior of F, and is denoted intF. Notice that there is only one type 
of interior. We define V^(F) := |intF|. The union of the remaining components is 
called the exterior of F, and is denoted by extF. A contour is external if it is not 
contained in the interior of another contour. 

Let F be a contour of some configuration a. Assume F is of type +. Consider 
the configuration o"[F], which coincides with on the support of F, and which 
equals on extF, —1 on intF. Using Proposition 12.11 it is easy to see that a[T] 
has a single contour, which is exactly F. This can be generalized to a family of 
external contours of the same type, as in the second part of the following lemma. 

Lemma 2.5. External contours have the following properties: 
1 ) External contours of an admissible family have the same type. 

^This property of R'' is called unicoherence. See |Kuj . vol. 2, Theorem 9 of Chapter 57.1, 
and Theorem 2 of Chapter 57.11. 
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2) Let {Ti, . . . , r^} be a family of external contours, all of the same type. Then 
{Ti, . . . , r.„} is admissible if and only if d{Ti, Tj) > I for all i ^ j . 

Proof. The first statement follows easily from Lemma I2.4I For the second, we 
can assume that the contours are of type +. If {Fi, . . . , F„} is admissible, then 
by construction the Fj are at distance at least I. Then, assume (i(Fj,Fj) > / for 
all i 7^ j. Consider the configuration cr[Fi, . . . , Fjj], which coincides with ar^ on 
the support of Fj, which equals +1 on f^jextFj and —1 on [JjintFj. Then the 
contours of cr[Fi, . . . , F„] are given by {Fi, . . . , F„}. □ 

2.2. Re-formulation of the Hamiltonian. Consider a finite volume A G C^'^ 
with the pure +-boundary condition +ac G Q^'^. Let a\ G Qa- We set a : = 
cta+ac The hamiltonian with boundary condition +ac is defined by 

HA{a) = //a(^a+aO = Yl ^^^^"^'^ ^-j) + '"^'^^) ' (2-2^) 

{i,i}nA7^0 jgA 

where u{(Ti) = —hai, /i G M. Since we work in a finite volume, we will from now 
on identify I* (a) with I* (a) HA and A^{a) with A=^((j) fl A. The following lemma 
shows how the hamiltonian can be written in such a way that spins in correct 
regions interact via the functions and are subject to an effective external field 
U#. 

Lemma 2.6. Define the potential U'^^ai) := u{ai) + Suppose 
a A is such that I* {a) H d]^A = . Then 

HAia) = Hj,iaj*aj,.) + J2{ E w*ia.,a,) + ^ U*i^^)) ■ (2.29) 

# {ij}nA#^0 ieA# 

Proof. The proof is a simple rearrangement of the terms. Consider a pair {i,j} 
appearing in Ha{ci). Since (i(A+, A~) > R we have a certain number of cases to 
consider: 1) {i,j} C A'^. In this case, write 

(t>ij{cri, Gj) = wtj{(Ti, aj) + (t)ij{ai, +) + 0ij(+, (Jj) . (2.30) 

The second term contributes to [/""^(crj), the third to U^{aj). 2) i G A+, j E I* . 
In this case the third term contributes to Hj*{ai*ai*<:). 3) z G A+, j G A'^; in this 
case, (pij{+,aj) = 0. The other cases are similar. Notice that the case i G A^, 
j G A'^ never occurs since points of d'^A can only be {6, +)-correct. □ 

2.3. Peierls Condition and Isoperimetric Constants. We take a closer look 
at the term Hj* . Remember that contours are maximal /^-connected components 
of /*. For each contour F, a[T] and aj*ai*c coincide on [I*]r. Since d(T,T') > I, 
we can decompose 

Hj,{ai,ai.c) = J2Hr{(T[T]) (2.31) 
r 

= E(liril + E^Hrw)' (2-32) 

r ier 



17 



where the sum is over contours of the configuration o (contained in A), and where 
the surface energy is defined as 

l|r||:= <^^.-(^[r]„a[r],). (2.33) 

{i,i}nr^0 

The central resuh of this section is the following. 

Proposition 2.2. The surface energy satisfies the Peierls condition, i.e. there 
exists p = p{6, u) > such that for all contour T, 

||r||>p|r|. (2.34) 

The constant p is independent of ^ and is called the Peierls constant. 
We will need two lemmas. The first is purely geometric. 

Lemma 2.7. For any finite set A C Z'^ and for all Rq G N, there exists Aq C A, 
called an RQ-dtpproxwrnnt of A, such that 

(1) A c [AoU , 

(2) d{x, y) > Ro for all x,y & Aq, x ^ y. 



The second lemma is a property of the Kac potential. In |BZ2j . this property was 
called "continuity" for obvious reasons. 

Lemma 2.8. Let a e ^, i e TJ^ , # G {±}. Define 

K(^;#):= 5^0.,(#,a,). (2.35) 

Then there exists C2 > such that for all x, y, d{x, y) < R, 

|V;(x;#)-K(?/;#)|<C2^^|^. (2.36) 
Proof. The difference K-(x; if) — V„{y] #) can be expressed as follows: 

X] (0xz(#,ctO - 0y/(#,(Tz)) - ^ (j)yi[j^^ai) 

ISBrIx) lf^Bii{x)r\BR{y) l&Bniy) 

l^Bniy) I^Br^x) 

The first and last sum can be estimated as follows: 

sup (2.37) 



1&Br{x) 



, , /2R+ l\d^^d(x,y) 
< do, ( sup J{t)) (— ^) (2.38) 
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Since we are considering the step function, sup^ J(t) = 2 The middle sum 
vanishes ^, which finishes the proof. □ 

Proof of Proposition By Lemma f2. 31 there exists in the 2i?- neighbourhood of 
each C^''' C F a point j G F that is (5-incorrect for o"[F]. Let A be the set of all 
such points. We have F C [^]«+2_r- Let Aq be any 4i?-approximant of A. We have 
A C [Ao]4iJ, i-e. F C [AoJi+eK- Each j G Aq is (5-incorrect for (t[F] i.e. satisfies 

\w%{T,a[T],)\>~5. (2.40) 
Since \wf,^{T,(r\^]k) \ = 20jfe(±, a[F]fe), 

V;[r](j;±) = 'PA±^^[^]k) > ^- (2.41) 
We bound the surface energy from below as follows: 



2 

jeAo keBR{j)nT i-.i^k 

= \ll E K[r](fc;a[F],)>i5^ 5^ K[r](A:;a[F],) (2.42) 

E ^-[r](^;^[ry (2.43) 

rf{fcj)<4l7R 

where C2 was defined in Lemma (2.81 Moreover we have, using ()2.36p . for each k 
of the sum, 

V,yr]{k- a[V]k) = K[r](j; ^[F]^) + (l4[r](A;; a[V]u) - 14[r](j; ^[F]^)) (2.44) 

~5 d{k,j) 5 6 _5 
- 2~'^~^- 2"'^4^-4- ^^-^^^ 

We have used the fundamental fact that the correctness of a point j does not 
depend on the value taken by the spin aj. This gives the lower bound 

||F|| > l|Ao|l|S^^(0)|^ > ^|5^^(0)||S,+6i?,(0)r'|F| > p|F| . (2.46) 

□ 



^Here we use for the second time the fact that we are considering the step function p.l4|l . 
Nevertheless, if J is an arbitrary iiT-Lipshitz function: 

Y \(l)xi{#,(Ti) - (f)yi{#,ai)\< Kc^j'^ d{jx,jy) 

l£BR(x)nBniv) leBR{x)nBri{y) 

< Kc,^^\Bn{x)\^^ . (2.39) 
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Since the Peierls constant is uniform in 7, we will be able to study the van der 
Waals limit at fixed /3. Proposition 12.21 allows to define, for = 1,2,..., the 
following numbers called isoperimetric constants: 

K{N) := inf |k > : 1/(1)'^ < K||r||, for all T, V{T) > iv} . (2.47) 

These constants will play a crucial role in the construction of the phase diagram 
and in the study of non-analyticity. Some of their properties are given in the 
following lemma. 

Lemma 2.9. The sequence K{N) is decreasing and there exists positive constants 
c_ , c+ such that 

c_7 < inf K{N) < sup K{N) < c+7 . (2.48) 

N 

As a consequence, the following limit exists 

K{oo) := lim K{N) . (2.49) 

N^oo 

Moreover, there exists for all e > a sequence (T]\r)j^>i, limjv_^oo V(Xn) = +00, 
such that for N large enough. 



(1 - e)K(cx))||r;v|| < V{Tm)— < (1 + e)A'(oo)||rjv|| • (2.50) 

Proof. K{N) is decreasing by definition. For the upper bound, use the Peierls 
condition and Lemma f2.1(JI hereafter: for all P, 

^ < ^ = -7 = c+7 • (2.51) 

||i II p|i I pi pu 

For the lower bound, we explicitly construct a large contour of cubic shape. 
Fix N and take M G N so that Am = [-M;+MY n Z'^, A^v G C^^\ \Am\ > 
2N. Consider the configuration a defined by cTj = —1 if z G Am, aj = +1 if 
i G A^. Clearly, I* (a) contains a single contour Tm (of type +). Using ()2.6|) . 
IITmII < ITmI < 2/|(9j*'AA/| = 2h'R\diAM\- Taking M large enough guarantees 
|Am| > V{Tm) > ||Aa/|. This gives, since |(9i+Am| = 2d\A 



M\ 



IIP II > pip+A I - = c_7\/ Pa/ . 2.52 

||Im|| 2 2z/_rt Ioj^Aa/I Sdu 

The existence of the sequence (PAr)Ar>i follows from the definition of K{N) and 
from the existence of the limit K{oo). □ 

Lemma 2.10. Let B G C'^^\ and let A he a maximal R-connected component of 
B'^. Then 

\B\>\d+A\>l\A\^. (2.53) 

Proof. Consider the edge boundary 6~^A := {e = {i,j) : i ^ A,j ^ A^}, where 
means that i,j are nearest neighbours. Decompose 6~^A = i?i U ■ ■ ■ U Ed, 
where Ea is the set of edges of 6~^A that are parallel to the coordinate axis a. 
Suppose e = {i,j), i & A, j G A^. Since A is maximal, Cj'"* C B. Moreover, 

Te ■■= {j,J + (j - i),J + 2{j - z), . . . ,i + (^ - l)(j -z)}cB. (2.54) 
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For all e, e' G E^, fl Te' = 0. So for all a, 

\dM\>\ U ^e| = $^|Te| = ^|i?.|. (2.55) 

Considering the inequality \S~^A\ < d maxo, | i^^ | and the standard isoperimetric 
inequality \S^A\ > 2d\A\~d^ finishes the proof. □ 

3. Restricted Phases 

Restricted phases intervene when a set of contours {F} is fixed (with a config- 
uration (Tr on each of them) and when we re-sum over all the configurations that 
have this same set of contours. The set of configurations having the same set of 
contours was completely characterized in Proposition 12. II We are thus naturally 
led to consider systems living in a volume A with a boundary condition ?7ac, with 
the constraint that each point i E [A]/j must be ^-correct. Our aim is to obtain 
a polymer representation for the partition function of such systems, and to show 
that the associated free energy behaves analytically at h = 0. As will be seen, 
the presence of the constraint will allow to treat the system in a way very similar 
to a high temperature expansion. The study of restricted phases we present was 
invented by Bovier and Zahradnfk in |BZ2j . At a few points our development 
differs slightly from theirs, so we expose all the details. 

A source of complication will be that the definition of polymers, as well as their 
weights, will depend on the boundary conditions specified outside A. Typically, 
the A we want to consider is the volume between a given set of contours and the 
boundary of a box. That is, the boundary condition is specified partly by the 
spins on the contours and partly by the boundary condition outside the box. To 
have an idea of the objects we are going to construct, see Figure El 

We will only treat the case -|-, the case — being similar by symmetry. Fix 
< 5 < 5 < 2~'^. Consider any finite set A G C*^'-*. First of all, we must consider 
boundary conditions of the following type: 

Definition 3.1. A boundary condition rj^c g VL\c is +-admissible if each i G [A]r 
is {5,+) -correct for the configuration +a%c. 

More intuitively, a +-admissible boundary condition means that when looked 
from any point i inside of A, there is a majority of spins -|-1 on the boundary. In 
our case (i.e. with the step function), this can be formulated as: for each i G [A]/j, 

\B],{^nB\<\\Bn{i)\^ (3.1) 
where the set B is defined by 

B = B{r]Ac) := {i e A' : (va^), = -1} . (3.2) 

In this sense, these boundary conditions are "good" ; there is hope in being able 
to control the -f-phase in the volume A. Notice that the boundary condition 
specified by a contour on its interior is always admissible. This is the reason why 
the parameter 6 was introduced in their definition. 



21 



We define the function tliat allows to realize tlie constraint obtained after Propo- 
sition ()2.1|) : consider a +-admissible boundary condition r]\c G Qa^. Let i G [A]^, 
o"A G r^A, and define 

J 1 if z is (5, +)-correct for aA^A- ,^ „^ 

li(<7A) := < . (3.3) 

I otrierwise . 

Tlien define 

1K) = V(aA):= n (3-4) 

Notice tliat 1(+a) = 1 since t/a^ is +-admissible. Tlie hamiltonian we use for tlie 
restricted system is the one obtained after the re-formulation of Lemma f2. 61 in a 
region of +-correct points. Set a := OATiA'^- The restricted partition function with 
boundary condition riA'^ is 

Z/(A;?7aO:= lMexp(-/? ^ w%{cy,,cy,) - ^Y.^^^^^^ ■ (3-5) 

We will show that Z^^ can be put in the form Z^^ = e'^^^^^Zr, where Zr is 
the partition function of a polymer model, having a normally convergent cluster 
expansion in 

iJ+ = {/i G C : Re/i > -|} . (3.6) 

The reason for logZ^"^ to behave analytically at /i = is that the presence of 
contours is suppressed by l(crA), and that on each spin cij = — 1 acts an effective 
magnetic field 

[/+(-!) = /i+^0,, = l + /i, (3.7) 
which is close to 1 when /i is in a neighbourhood of /i = 0. 

3.1. Representation with Polymers. The influence of a boundary condition 
can always be interpreted as a magnetic field acting on sites near the boundary. 
We thus rearrange the terms of the hamiltonian as follows: 

Y wtj{^^,^j) + E (^^(^^) + E (^aO.)) • (3.8) 

{ij}CA «GA iGA<= 

By defining an new effective non-homogeneous magnetic field 

fitia,) := U+{a,) + h+Y (^a^)^) ' (3-9) 

we can extract a volume term from Z^"*" and get Zj,'^ = e^'^^^^Zr, where 

Z/:= J2 lMexp(-/3 Y n,t^icT,,a,)-l3Yf^ti^i))^ (3-10) 

ctaGSIa {ijjcA ieA 

Notice that the field ^^(cTj) becomes independent of ?7ac when d{i, A'^) > R. Since 
w^j{ai,aj) = if CTj = +1 or ctj = +1 and /x^(+l) = 0, we need only consider 
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points i with cTj = —1, which will be identified with the vertices of a graph. Each 
vertex of this graph will then get a factor e~^^» When h G 

Re^+(-l) = 1 + 2Re /i + ^ w±{-, (r/AOi) > 1 " 2| - ^ > ^ • (3.11) 
We used the fact that 6 < 2~'^. 

The formulation of Zr in terms of polymers will be a three step procedure. We 
first express as a sum over graphs, satisfying a certain constraint inherited 
from l(crA). Then, we associate to each graph a spanning tree and re-sum over all 
graphs having the same spanning tree. We will see that the weights of the trees 
obtained have good decreasing properties. Finally, the constraint is expanded, 
yielding sets on which the constraint is violated. These sets are linked with trees. 
After a second partial re-summation, this yields a sum over polymers, which are 
nothing but particular graphs with vertices living on Z*^ and whose edges are of 
length at most R. 

3.1.1. A sum over graphs. Let Q\ be the family of simple non-oriented graphs 
G = {V, E) where V d each edge e = {i, j} e E has d{i, j) < R. For e = {i, j}, 
set w'^ := wfj{—,—). Notice that = —2(j)ij < 0. Define also fif := 
Expanding the product over edges leads to the following expression 

Zr=j2 n (^'^""^ - 1) n ^"""^ ' (3.12) 

Ge^A eeS(G) i£V{G) 

where 1{V) := l{aA{V)), and a^iV) E is defined by aA{V)i = -lifiEV, +1 
otherwise. With this formulation in terms of graphs, the constraint 1{V{G)) = 1 
is satisfied if and only if for each i G [A]/?, 

e={«j} 

jeV{G)uB 

Moreover, the fact that the boundary condition r]\c is +-admissible reduces to 

e={«j} 

3.1.2. A sum over trees. Suppose we are given an algorithm that assigns to each 
connected graph Go a deterministic spanning tree T{Go), in a translation invari- 
ant way (that is if Gq is obtained from Go by translation then T'(Go) is obtained 
from T'(Go) by the same translation). To be precise, we consider the Penrose 
algorithm considered in Chapter 3 of jPj ^. We apply the Penrose algorithm to 
each component of each graph G appearing in the partition function. Let T\ G Q\ 
denote the set of all forests. We have 

Zr= Y.l{V{T))llu+{t), (3.15) 

tsTa ter 



The Penrose algorithm requires the choice of an origin among the vertices of the graph. We 
choose this origin as the smallest vertex of the graph with respect to the lexicographical order. 
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where the product is over trees of T, and the weight of each tree is defined by 



a;- t : = 



E n (^"''"' - 1) n ^-''^ ■ (3-16) 



GgSa: e<^E{G) i&V{G) 
T{G)=t 



Notice that isolated sites {i} C A are also considered as trees. In this case, 
uj^{{i}) = e~l^^i . The following lemma shows how the re-formulation in terms of 
trees allows to take advantage of the constraint. 

Lemma 3.1. Let T & be a forest such that 1{V{T)) = 1. Then for each tree 
teT, 

ii^+(t)b, < n n (3.17) 

eeE(t) i£V(t) 

Proof. For each t G T, let E*{t) denote the set of edges of the maximal d- 
connected graph of {G G Qa : T{G) = t} (see jPj)- We can express the weight 
as follows: 

a;+(t)= H (e-^^^-l) n e'^^' Yl n (3-18) 

eeEit) iey(t) GgSa; e£E{G)\E{t) 

T{G)=t 

= Yl (e-^"'^ - 1) n n ^'^^^ ■ (3.19) 

e6_E(t) iey(t) e£E*(t)\E{t) 

Since l(y(T)) = 1, the constraint ()3.13j) is satisfied, and the last product can be 
bounded by: 

II e('^^^^< n n ^^'"'*^' (3-20) 

ee£;*(t)\E(t) iey(t) e={ij} 

iey(t) 



iev'(t) e={i,j} ievit) 

This gives the result, since Re/i^ > | by (|3.1H) . and 6 < 2^'^ < j. □ 
Notice that to obtain ()3.2H) . we only needed that the bound 

\^t\ < ^ V? G V{t) (3.22) 

e={i,j} 
j^V(t) 

be satisfied. This is weaker than (j3.13|) and clearly l(y(T)) = 1 only if (j3.22j) 
is satisfied for all t G T. In the sequel we can thus assume that the trees we 
consider always satisfy (j3.22|) . independently of each other. So the bound (j3.17|) 
can always be used. A direct consequence of the last lemma is the following result 
which shows that trees and their weights satisfy the main condition ensuring 
convergence of cluster expansions. 
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Corollary 3.1. Let < c < e > 0. There exists 7o > and (3i = Pi{e) such 
that for all 7 G (0,7o), (3 > jSi, the following hound holds: 

J2 lk+(t)|k+e^l^(*)l <e. (3.23) 

t:y(t)90 

Proof. Using Lemma 03 

lk+(t)||^,^e^l^(*)l < Yl {e-"""^-!) Yl e-'s". (3.24) 
eeEit) ieV{t) 

When t is a single isolated point (the origin), then we have a factor e~^^. When 
V{t) 3 0, E(t) 7^ 0, we define the generation of t, gen(t), as the number of edges 
of the longest self avoiding path in t starting at the origin. The sum in ()3.23p is 
bounded by 

+ E E n (^"^"^^ - 1) n (3.25) 

<7>1 t:y(t)90 eG£;(t) i€V{t) 
gcn(t)=g 

<e--sf' + J2e-^^''' Yl n (e^^'-^-l) H ^""^^ (^.26) 

9>1 t:y(t)90 eGS{t) jeV(t) 

gen(t)=g 

< e-^'^ + ^e-ii'^%g, (3.27) 
where we defined {Vi{t) is the set of leaves of the tree t): 

E n (^"""' - 1) n n ^-^^ (3-28) 

t:V(t)50 eeE{t) ieV{t)\Vi(t) iey;(t) 

gcn(t)=g 

Before going further, we define 

70 := sup {7 > : 20^7"^ sup J{s) < ^} . (3.29) 

Since e~^'^^ — 1 < /?|i«^|e^'"''^' and {w^l = 20jj we can bound, when 7 < 70, 

(e"^"=^ - l)e-3^^ < /3e-^^ ^ < 2(3e-^^^ = /3C(/3) . (3.30) 

We are going to show that a^+i < ag for all (7 > 1. Clearly, a tree t of generation 
g + 1 can be obtained from a sub-tree t' C t of generation g by attaching edges 
to leaves of t'. Let x be a leaf of t'. The sum over all possible edges (if any) 
attached at x is bounded by 

1 + E H ^ ■ ■ ■ E n (<=""•"•• - 1)-"*" s 1 + E B^i^mf = . 

k>l eiSx e^Bx i=l k>l 

Assuming (3 is large enough so that ({(3) < ^, the weight of the leaf x changes 

into e~^^e^'^^^^ < e~^^, which is exactly what appears in Og. This shows that 
CKg+i < C(g. We then have a^+i < ctg < ■ ■ ■ < ai. Like we just did, it is easy to 
see that cti < e~32^. This proves the result. □ 
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3.1.3. A sum over polymers. After the partial re- summation over the graphs hav- 
ing the same spanning tree, the constraint 1{V{T)) in ()3.15p still depends on the 
relative positions of the trees. This "multi-body interaction" can be worked out 
by expanding 

Hv{T))= n uv{T))= i[{i + imm= E U^^inT)), 

«g[a]r je[A]R MclA] RieM 

where l^(y(r)) := li{V{T)) - 1. This yields 

^r=J2 H (ni^(^(^)))(n'^^W)- (3-31) 

T&Ta Mc[A-]r ieM t6T 

Consider a pair {T,M) in (j3.3ip . Let i G M. The function l^(y(T)) is non-zero 
only when i is not {6, +)-correct; it depends on the presence of trees of T in the 
i?-neighbourhood of i and eventually on the points of B{rj\c) if Buij) fl A'^ 7^ 0. 
To make these dependencies only local, we are going to link the /^-neighbourhood 
of points of M with the trees of T. 

Consider the graph = N{M) defined as follows: the vertices of A^ are given by 

ViN) := \J Br{i) . (3.32) 

iGA/ 

Then, A has an edge between x and y if and only if (x, y) is a pair of nearest 
neighbours of the same box Bjiii) for some i G M. The graph A decomposes nat- 
urally into connected components (in the sense of graph theory) Ai, A2, . . . , Nk- 
Some of these components can intersect A'^. 

We then link trees tj G T with components A^ G A^. To this end, we define 
an abstract graph G: to each tree tj G T, associate an abstract vertex wi and 
to each component Nj an abstract vertex Zj. The edges of G are defined as fol- 
lows: G has only edges between vertices Wi and Zj, and this occurs if and only 
if V{ti) n V{Nj) ^ 0. Consider a connected component of G*, whose vertices 
{wi^ ,...,Wi^,Zj^,..., Zj^ } correspond to a set P/ = {tj^ Nj^, Nj^}. We 
change P[ into a set Pi, using the following decimation procedure: if P/ = {tjj} 
is a single tree then Pi := P/. Otherwise, 

1) delete from P/ all trees tj^ that have no edges, 

2) for all tree tj^ containing at least one edge, delete all edges e G E{ti^) whose 
both end-points lie in the same component Nj^. 

The resulting set is of the form Pi = {t^-^, . . . , t^^, A^^, . . . , Nj^}, where each tree 
ts. is a sub-tree of one of the trees {tj^, . . . ,tjj}. P; is called a polymer. The 
decimation procedure P/ =^ P/ is depicted on Figure El 

The body of Pi is B{Pi) := V{NjJ U ■ ■ ■ U V{Nj^). The legs of Pi, C{Pi), are the 
trees {t^^, . . . ,tsj. 

A polymer can have no body (in which case it is a tree of 7a), or no legs (in which 
case it is a single component AjJ. We define the support V{P) as the total set 
of sites: 

V{P):= U V{t)u[jV{N,). (3.33) 

te£(P) * 
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t 




Figure 4: The decimation procedure P/ Pi. The hatched polygons represent the 
body B{Pi) and the legs are the trees {tg^ , t^^ , tsg, t^^, ts^}. Each tg^ is a sub-tree of 
some ti. 



Often we denote V{P) also by P. Two polymers are compatible if and only if 
V{Pi) n V{P2) = 0, denoted Pi ~ P2. We have thus associated to each pair 
(r, M) a family of pairwise compatible polymers {P} := (p{T, M). The set of all 
possible polymers constructed in this way is denoted Vj^lrj/^c). The representation 
of Zr in terms of polymers is then 

J2 n '^^(^) ' (3.34) 

{P}CV+{riAc) Pe{P} 
compat. 

where the weight is defined by 

oo^p) E ( n ( n ^^(^)) (3.35) 

(T,A/): i&M teT 
ip{T,M)=P 

We should have in mind that uj^{P) depends on the position of P inside the 
volume A, via the boundary condition r]\c: more precisely if B{P) fl A*^ 7^ or 
if there exists a leg t G C{P) such that d{t,A'^) < R. Therefore, we define the 
family of free polymers whose weight depends only on the intrinsic structure of 
P, and not on the boundary condition. The family is translation invariant, 
as well as the weight of each of its polymers. To any finite family V, we associate 
the partition function 

MV) := E n ^^(^) ' (3-36) 

{P}CV P6{P} 
compat. 

where the product equals 1 when {P} = 0. For instance, we have obtained 

Z/(A; r^Ac) = e^'^I^IZ,,(P+(77Ae)) . (3.37) 

Everything we have done until now can be done for a — admissible boundary 
condition tac, yielding a family of polymers 'Pa(ta<:), with weights lu^{P). In 
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this case, sites get a factor e'^'^^ . In particular, if we consider the spin-flipped 
boundary condition — r/A^ defined by {—r]/<^c)i := — (t^aO*) which is — admissible, 
we have when h is purely imaginary ^, 



MVUVA^)) = ZriVXhVA^)) . (3.38) 
This symmetry holds only because we consider free polymers. 

3.2. Analyticity of the Restricted Phases. Define the restricted pressures by 

p% := hm --i^logZ,±(A;±AO , (3.39) 
' A/'Z'* p\A\ 

where the thermodynamic limit is taken along a sequence of cubes. A result of 
the present section is that the restricted pressures, unlike the total pressure p^, 
behave analytically at /i = 0. 

We study the weight a;"^(P) {uj^{P) is similar by symmetry). The point is that 
we linked trees with the i?-neighbourhood of points of the set M, and we must 
now see that this thickening does not destroy, from the point of view of entropy, 
the uniformity we have been able to obtain with respect to the scaling parameter 
7. Moreover, the body of polymers can intersect A'^. At this point we will see 
that 5 — (5 > is crucial. 

Lemma 3.2. There exists (32 O'nd tq > such that for all f3 > one? for all 
7 G (0,7o), the following holds: each polymer P G V^i^rj/^c) satisfies 

lk+(P)||H+<e--«^l^(^)l n ie-^""^-!) H e-i^^ (3.40) 

ee£(P) iG£(P) 

Proof. Remember that the bound (j3.17|l holds for each tree under consideration. 
If B{P) = 0, then P is a tree and the result follows from Lemma [3. II Otherwise, 



11-^^)11^. < E ( n i^'iinm) n ( n (^"''"' - 1) n ^"^') • 

(T,A/): i&M t&T eG-B(t) if^V{t) 

ip{T,M)=P 

Consider a pair (T, M) such that if{T, M) = P. Let zq G M, and assume 
l^^(y(T)) 7^ 0. This implies, with regard to ()3.13|) . 

E l^e'-lx^- (3.41) 

V 

But, according to ()3.14|) . we have 



e={io,j} 
jeV{T)UB 



e={io,j} 



This implies the crucial lower bound 

J2 \wt\>S-6>0. (3.43) 



e={io,j} 



Here, z denotes the complex conjugate of z. 
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Since \w^\ = 2(pij < 20^7*^ sup j, J{s), we can find a constant C3 such that 

\V{T) n BUzo)\ >{6- ~S)cs\Bn{io)\ . (3.44) 

In this sense, the forests that contribute to u}^{P) accumulate in the neighbour- 
hood of each point zq € M. See Figure El Let Mq be any 2i?-approximant of M. 
Then we have \B{P)\ < |Mo||S3r(0)| and so 

\V{T)nBiP)\ > J2 \V{T)nBnito)\ > (S - ~6)c,\BiP)\ (3.45) 

where C4 is a constant. Now, each i G V{T) gets a factor e~^^ = e~'^^^. One 
factor e~i2^ contributes to extract a term decreasing exponentially fast with the 
size of B{P), using ^J^: 

A second factor e~A^ contributes to the weight of the legs. Extracting this 
contribution gives 

-Q -Q g-^/?^ (3 47) 

eeC(P) i(iC{P) 

The last factor e~^^ is used to re-sum over all the possible configurations of T 
inside the body B{P) (see Figure 13)), that is over all forests T', V{T') C B{P), 
where each tree t' G T' gets a weight bounded by 

uJo{t'):= n {e-^'"^-l) n (3.48) 

The remaining sum is thus bounded by 

J2 n = 0o(S(P)) . (3.49) 

T':V{T')CB{P) t'eT' 




Figure 5: The re-summation of Lemma f3. 21 We emphasized the fact that the forest 
T must have many points in B{P) n A, as was shown in (|3.45|l . 
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This partition function can be studied with a convergent cluster expansion. Pro- 
ceeding as we did in Corollary 13.11 we can take (] sufficiently large so that the 
weight ujoit') satisfies ()3.23|) . We can then guarantee that 

|logeo(S(P))|<|i3(P)|. (3.50) 

The sum over all possible sets M such that N{M) has a set of vertices given by 
B{P) is bounded by 2^'^^^'^^. Altogether these bounds give 

g-i{5-5)c4/3|B(P)|2|B(P)|g|B(P)| = g-ro/3|B(P)| ^ 

which finishes the proof. □ 

We now give the consequence of this lemma, namely that polymers satisfy the 
main criterion needed for having a convergent cluster expansion. 

Corollary 3.2. Let < c < min(^, e > 0. There exists (3^ = /93(e), such 
that for all P > Ps and for all 7 G (0, 70) , the following holds: 

Y: |k^(P)||H,e^l^(^)l<e. (3.51) 

P:y(P)90 

Proof. Lemma f3.2l allows to bound 

||^+(P)b+ < ( n ^o(iV)) ( n ^o(^)) = ^o(^) ' (3-52) 
NeP t6£(P) 

where the weight of each component of the body is ujo{N) := e~'^°'^l^(^)l and 
the weight of each leg t was defined in ()3.48|) . Fix e > small. It is easy to show 
that when P is large enough, 

J2 ^o(iV)e('=+^)l^(^)l<|e, (3.53) 

N:V{N)30 

and, proceeding like in Corollarv 13.11 

J2 a;o(t)e(^+^)l^(*)l < |e. (3.54) 

t:y(t)30 

Let n{P) denote the number of objects (components N and trees t) contained in 
P. That is, if P = {ti, . . . , ti, Ai, . . . , Nk}, then n(P) = L + K. We will show 
by induction on A = 1, 2, ... ^ that 

A;v:= Yl ^0(^)6^=1^(^)1 <e, (3.55) 

P:V(P)30 
n{P)<N 

which will finish the proof. If A = 1 then P can be either a single component A 
or a tree t. The bound then follows from (|3.53p and (|3.54j) . Suppose P is large and 
that the bound holds for A^. If P satisfies V{P) 3 0, n{P) < A + 1, we choose an 
object of P that contains the origin (which can be a tree tg or a component Aq), 
and decompose P as follows: either P = {Aq} U {Pi, . . . , P^} with V{Nq) 3 0, 
V{P,) n VjN o) ^ 0, n(P,) < A, P, ~ P,- for z ^ j, or P = {to} U {Pi,...,Pfc} 

^We thank Daniel Ueltschi for suggesting this method, which wih be used again later in the 
study of the phase diagram. 
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with V{to) 3 0, and V{P,) n V{to) ^ 0, n{Pi) < iV, ~ Pj for i j. In the 
first case, we have, using the induction hypothesis and ()3.53|) . 

Y: ^o(A^o)e^l^(^«^IE^( E -o(P)e^l^(^)l)' (3.56) 

No:V{No)50 k>0 ' P:V{P)nV{No)j^(/> 

n{P)<N 

< c.o(iVo)e<^l^(^°)IE^(l^(^o)|A^)^ (3.57) 

No:V{No)30 k>0 

< Y ^o(iVo)e'^l^(^'')le^l^(^°)l<|6. (3.58) 

No:ViNo)90 

In the second case the same computation yields, using ()3.54|1 . 

Y: -o(to)e^l^(-)l$:i( Y -o(F)e^l^(^)l)' 

to:V(to)30 A:>0 ' P:V{P)nVito)^9 

n(P)<N 

< Y ^o(to)e^l^(''')le^l^(*")l < |e. (3.59) 

This shows that Aat+i < e and finishes the proof. □ 

We now state the main result concerning restricted phases and their analyticity 
properties, again only for the case # = +. We refer to Appendix ^ for nota- 
tions. Here polymers play the role of animals. Clusters of polymers associated to 
'P][{'r]\c) are denoted P G Vj^i^ij^c). By Lemma Qj. 11 (j3.5ip implies 

supV||u;+(P)||^^^ <supV|a;o(P)| <?7(e), (3.60) 

P^x P^x 

where the weights uj^{P) and uJo{P) are defined like in ()B.3|) . Since e can be made 
arbitrarily small by taking (3 large enough, we will replace rjle) by a function er{(3), 
where the subscript r is to indicate that this function concerns the restricted 
phase. We define if+ := {Reh > — ^} C 

Theorem 3.1. Let /3 be large enough, 7 G (0,70). Let A G C-'^ and r]\c be a +- 
admissible boundary condition. Then Zr{VX{f]h'')) has a cluster expansion that 
converges normally in given by 

logZ,,(P+(r;AO)= Y ^''(^)- (3.61) 

The maps h 1— > \og Zr{V^{r}\c)) , h h-* p^^{h) are analytic in Moreover there 
exists a function er{(3), lim^yoo er(/3) = 0, such that 

II logZ,,(P+(77Ac))||^^ < eriPMl , Y < 'riP) , (3.62) 

||^iogZ,(P+(r/AO)L, <e.(/?)|A|. (3-63) 
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The proof of the theorem follows easily from Lemma IB.ll Analyticity follows 
from the fact that the convergence is normal on H^. The bound on the first 
derivative is obtained by using the Cauchy formula: any disc of radius centered 
at 2; e if+ is contained in if^.. This also implies the existence of a constant > 
such that for all integer k > 2, 

logZ/(A; ?7aO < C^kl , (0)1 < C^kl . (3.64) 

h=0 

4. The Phase Diagram 
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Throughout this section and until the end of the paper we assume 7 G (0, 70) is 
fixed, where 70 was given in ()H.29j) . To start with, consider the partition function 

Z+(A) := J2 e-^^^("^+^^) , (4.1) 

where 

nt := {^A e : d(r(aA+AO, A^) > /} . (4.2) 

For each a a G Qj^, the decomposition of /*(o"a+ac) into connected components 
yields an admissible family {F}, such that F C A and d(r,A'^) > I for each 
F G {F}. Then, A is decomposed into A = {F} U A"*" U A~, where A* are the 
points of A\{F} that are (5, #)-correct for the configuration a"A+A'=- 
In fj4.1|) . we re-sum over the configurations a\+ (resp. cr\+) on A"*" (resp. A") that 
yield the same set of contours {F}. In Proposition 12 . 1 1 we characterized explicitely 
the constraints satisfied by the configurations aA±: each point i G [A"*"]/; must be 
{6, -|-)-correct for the configuration o"a+ -I-a'= ^"{r}, where a^r} is the configuration 
specified by the contours on the union of their supports. Similarly, each point 
i G [A^]r must be (5, — )-correct for the configuration a\-o-{r}- Using the re- 
formulation of the hamiltonian given in Lemma f2. 61 we get: 

Z+(A)= J2 ( n P(r))^/(A+;+Aca{r})^r(A-;a|r}), (4.3) 
{r}cA re{r} 

where the sum is over admissible families of contours, and 

p(F) := e-^^r(-[ri) . (4.4) 

Notice that when {F} = 0, then A = A''" and the summand of ()4.3p equals 
Zr~^{A; +ac)- Since they are subject to boundary conditions that depend on the 
family of contours {F}, the restricted phases induce an interaction among the 
contours. Nevertheless, the boundary conditions imposed by the contours and 
+\c on A"^ and A~ are admissible (in the sense of Definition 13. Ij) . This implies 
that the results of Section 01 can be used for the restricted partition functions 
appearing in (|4.3|) . 



Since we need to represent the partition function with objects whose compati- 
bility is purely geometrical, we need to proceed by induction, and consider systems 
living in the interior of external contours. Therefore, we must study functions 
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similar as ()4.3|) . with an arbitrary +-admissible boundary condition rjAc. We thus 
define 

0+(A;r/Ac):= ^ ( J] p(r))z/(A+; r/Aca{r})^r(A'; ^{r}) , (4.5) 
{r}cA re{r} 

Contours always lie at least at distance / from A'^. The external contours contours 
of {r} can be subject to particular constraints (as will appear, for example, in 
Section|S)), but we omit it in the notation. Notice that the empty family {F} = 0, 
the summand corresponds to a pure restricted phase Zr~^{A; r]\c). 

The aim, in the study of 0"'"(A; ?7ac), is to extract from (|4.5|) a global contribu- 
tion of the restricted phase. In the Ising model, the same operation amounts to 
extract the trivial term e^'*'^'. Here we extract Zj,'''(A, t^a^^) = ^^'^''^'-^^(Pa (%''))' 
and our aim is to reach the representation ()4.17p . The deviations from the re- 
stricted phase will be modelized by chains, i.e. contours linked by clusters of 
polymers (polymers describe the restricted phase). In Section H^Tl we expose this 
linking procedure. In Section we show how to handle the entropy of chains, 
preserving the uniformity in the scaling parameter 7. In Section we study the 
weights of chains and their dependence on the magnetic field near Re = 0, i.e. at 
coexistence. In Section lOl we study pure phases, i.e. {Re/i > 0} and {Re/i < 0}. 



4.1. The Linking Procedure. We first express Q^{K\ri\c) as a sum over ex- 
ternal contours. By Lemma (2.51 each external contour is of type -|-. Let {F} be 
a family of external contours. Then, A is decomposed into 

A = extA{F} U {F} U |J intF , 

re{r} 

where extA{F} := A n Hreir} G^tF. For each family of admissible external con- 
tours {F}, we re-sum over the configurations whose external contours are given 
exactly by {F}. This induces, for all F, a partition function 6~(intF; +o"r), which 
can be expressed as in fl4.5p . On extA{F}, we get a restricted partition function 
Zr'''(extA{F}; r^A'^o'ir}). We thus have 

e+(A;77AO = 

Z/(A;?7aO+ ^r+(extA{F};r/Aca{r})nP(r)0"(intF;ar), (4.6) 
{r}cA r 

ext. 

where the sum is over non-empty families of external contours. Consider the 
configuration —err obtained by spin-fiipping ar, i.e. (— crr)j := — (c"r)j for all z G F. 
We introduce the functions Z,.'''(intF; — cxr) and 6"'~(intF; —err) and consider, for 
a while, 

Z^+(extA{F}; r^A'^air}) Hr ^r^(intF; -ar) 
Zr+(A; riKc) 

Using the polymer representation of Section IHl we consider the family of poly- 
mers := Pg^^ |r|('7A''^{r}) associated to Zr~^(extA{F}; 77A'=cr{r}), the families 
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^intr ^intr(^^r) associated to each of the ^^"''(intr; — ar), as well as the family 
VX '■= T^XijIh-'') associated to Zr'^(A; r]i^c). Since the expansions of these functions 
are absolutely convergent, we can rearrange the terms. The volume contributions 
from extA{r} and IJp intF cancel, and we get 

where we used the abbreviation 

5^±^+(P)= u^P)- Yl ^^(^)- (4-8) 

p p&v+t P&rl 

d(P,{r})<R d{P,{r})<R 

pncxtA{r}7^0 

The sign ± in front of uj~^{P) is chosen in function of the sum to which P belongs. 
Define A+(P) := e^'^^^^) - 1 and expand 

eEp±^-(P) = J|(l + A+(P))= Y nA+(A). (4.9) 

P {Pl,-,Pn} * = 1 

The function depends only on the structure of F, and is given by 

E^= Y ^^(P)- E ^^(^)' (4.10) 

d(P,r)<R Pnextr=0 
d{p,r)<R 

where P"*" denotes the family of clusters associated to free polymers of type +. 
Notice that is analytic in H^. Since |[r]^| < 3'^|r| we have, if /? is large 
enough (see Theorem 13. 1|) 

ll^r lk+ < -|r| , ll-TT^rllff <-|r|- (4.11) 

If we define the weight (we denote +ar = or) 
with pi{T) := p(r)e-^''l^le^?, we have 

ext. 

We can then repeat the same procedure of summing inside external contours 
of 0'''(intr; — o"r), etc. This procedure continues until we reach contours whose 
interior can't contain any contour. At the end, 

H^-i+E5:(n^^(»)(n-^(r)). (4.13) 

^ ' ' {r}cA {P} p r 
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where the sum {F} contains contours of type +, and each cluster P lies at distance 
at most R from one or several contours of {F}. For this reason, the weight of 
some polymers can depend on the configuration ar of the contours F that lie in 
their neighbourhood (or on r]\c). 

We get rid of these dependencies by linking polymers to contours. Like we did 
in Section 121 (when linking trees with components of the graph N), we associate 
to each pair ({F}, {P}) an abstract graph G as follows: each contour Fj G {F} 
is represented by an abstract vertex zj, each cluster P^ G {P} is represented by 
an abstract vertex Wk- This defines V{G). Then, we put an edge between zj and 
Wk if and only if d(Tj, Pk) < R. We also put an edge between Wk^ and Wk^ if and 
only if l^(Aj n V(Aj ^ 0. 

Each connected component of G, with vertices, say, {zj^^, . . . , Zji,Wki, ■ ■ ■ ,Wki}, 
represents a subset of {F} U {P} given hj X = {Fj^, . . . , Fj^ P^^, . . . , Pfc^}. X is 
called a chain of contours, or simply a chain. We denote by {X} the family of 
chains associated to the pair ({F}, {P}). The chains of {X} are of type +, and 
pairwise compatible by definition. The support of X, also written X, denotes the 
union IJrex -'^ Upex-^- Notice that if two chains X,X' are not compatible, 
then b{X) n b{X') ^ 0, where 

b{X):= [jnu U P. (4.14) 

rex pgx 

The weight of a chain is defined by 

u^'-iX) := ( n ^^(P)) ( n ^^(n) , (4.15) 
Pex rex 

and depends only on the intrinsic structure of the chain X (except, maybe, if 
d{X,A'^) < R). The final representation of the partition function is thus 

e+(A;r/Ac) = Z/(A;7;Ac)^ J] cu+(X) (4.16) 

{X} xe{x} 

= Z,+(A;?7Ac)S+(A;?7Ac). (4.17) 

In ()4.16p . the product is defined to be equal to 1 when {X} = 0. This last 
expression nicely expresses the fact that chains of contours describe deviations 
from a restricted phase. For the restricted phase, there corresponds a family 
Vj^{r]\c) associated to Zr~^{A] rjt^c). Similarly, there corresponds a family of chains 
'^A (^A=) associated to S+(A; r]]yc). The partition function can be written in terms 
of these families as 

e+(A; t^aO = e^^l^lZ,(P+(r7A0)2(A'+(r;A0) • (4-18) 

By definition, 'E{X'l^{rji^c)) := 1 when Xj^rj^c) = 0. Everything that was done 
until now can be applied also to the case where ri\c is — admissible, yielding 
chains of type — . 
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4.2. The Entropy of Chains. Before starting the analysis of the weights, we 
show how a priori bounds on the weights X^{P) and cu"*'(r) allow to handle the 
summation of weights of chains. In this section we assume that |A^(P)| < Ao(-P), 
|a;+(r)| <po(r),i.e. 

< ( n ^o(p)) ( n po(n) = ^o(^) . (4.19) 

Pex rex 

Convention: Now and in the sequel we will always use a subscript "0" in the 
weight of an object to specify that it depends only on the geometric structure of 
the object (as we did in ()3.52j) . Section IT^ . That is, such weights will always be 
translation invariant. When a weight is defined for an object, we use the same 
letter for the weight of the clusters of such objects (see Appendix IB|i. 

The proof of the following lemma is essentially the same as the one of Corollary 
EH We use the notations |F| := | [Jp^pV{P)\, \X\ := ^p^^ |r| + J2pex 1^1- 

Lemma 4.1. Let c > 0, e > 0, and assume the weights Ao(-P), Po(X) satisfy the 
bounds 

J2 Ao(F)e(^+^(^^+^)l^l < I , Y Po(r)e(^+^)ir]'l < . (4.20) 

Then the weight ujq{X) satisfies the condition ()B.4|) of Lemma \B . 1\ Namely, 

Y cuo(X)e'=l''(^)l < e. (4.21) 

X:6(X)90 

Proof. For a chain X = {Fi, . . . , Fx, Pi, . . . , Pm}, let n{X) := L + M denote the 
number of objects composing X (a cluster is considered as a single object). We 
show by induction on = 1, 2, . . . that 

^n:= J2 ^o(X)e^l''WI < 6 (4.22) 

X:fe(X)90 
n{X)<N 

If n{X) = 1 then X contains a single object, i.e. a contour. Then .^i < e follows 
from ()4.20|) . So suppose ()4.22|) holds for N, and consider ^n+i] this sum can be 
bounded by a sum in which each chain X is decomposed into [Fq]; 3 0, X 3 Fq, 
or into Po 3 0, X 3 Pq. This means: 

1) in the first case, X decomposes into A = {Fo}U{Xi, . . . , Xi^} ^° with [Fq]; 3 0, 
d{Xi, Fo) < R, n{Xi) < A for aU z = 1, . . . , n Xj = for aU i ^ j. The 



The chains Xi are obtained as follows: consider the abstract connected graph G associated 
to the chain X. Then, remove all the edges of G that are adjacent to the vertex zq represent- 
ing To and zq itself, and consider the decomposition of the remaining graph into connected 
components. These components are exactly the representatives oi Xi, . . . , Xk- 
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contribution to ^tv+i is thus bounded by 

E /^o(ro)e<^l[^»''IE^n E -o(^.)e^l^(-^^l (4.23) 

ro:[ro]i90 K>0 ■ i=l Xr.d(X„Vo)<R 

n{X,)<N 

< E Po(ro)e^i[^«i'iE^(l[ro]^ie^)"< E Po(ro)e(-^)i[^oi'i < I , 

ro:[ro],90 K>o ■ ro:[ro],30 

where we used the induction hypothesis < 

2) in the second case, X = {Pq} U {Xi, . . . , Xk} with Pq 3 0, d{Xi, Pq) < R, 
n{Xi) < N for all i = 1, . . . , K, Xi n Xj = for all i ^ j. A chain Xi of this 
decomposition can be of two types: i) there exists a cluster P & Xi such that 
P n Pq 7^ 0. Then the contribution from these chains is at most 

|Po| E ^o(^^)e'=l''(^')l = iPoie^v < \Po\e. (4.24) 

X,:b{X,)BO 
n{X,)<N 

ii) there exists V G X^, V fl {[Po]_R}i 7^ 0, where the thickening {■}; was defined 
in (ITTH|) . Notice that the set {[Po\r]i e C^') contains at most 2'^|Po| cubes C^'). 
Since contours are composed of cubes C-'-*, the contribution from these chains 
can be bounded by 

2lPo|eiv<2"e|Po|. (4.25) 
We can then proceed like in (j4.23j) . and get a contribution to C,n+i bounded by 

E Ao(Po)e^l^°le^(''+^)l^°l < | (4.26) 

Altogether, this shows that ^n+i < e. □ 

4.3. Domains of Analyticity. In this section we consider the dependence of the 
weights uj'^{X) on the magnetic field e C, in a neighbourhood of {Reh = 0}. 
For obvious reasons, the domain in which uj~^{X) can be shown to be analytic 
depends on the contour T E X that has the largest interior. Everything we say 
in this section holds for chains of both types, but for the sake of simplicity, the 
statements will be given only for chains of type +. 

The domains of analyticity depend on the isoperimetric constants K{N) defined 
in (j2.47j) . Consider the reals 

R{N) := , (4.27) 

2K{N)N-^ 

where 9 G (0, 1) will play an important role later in the study of the derivatives. 
We know from Lemma f2. 91 that R{N)Nd is increasing and that 
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Since we want the domains of analyticity to be decreasing with the size of the 
contours, we define 

R*{N) := min {R{N') : 1 < N' < N} . (4.29) 

The sequences R*{N) and R{N) have the same asymptotic behaviour, as the 
following lemma shows. 

Lemma 4.2. 

lim R*(N)N^ = (4.30) 

N^oo ^ ' 2K{oo) ^ ' 

Proof. First notice that there exists an unbounded increasing sequence Ni, N2, . . . , 
such that R*{Ni) = R{Ni). This is a direct consequence of the bounds 

R*iN)<RiN)< \ . (4.31) 

Since R(N)Nd increases, it is sufficient to show that R*(N)N^ is increasing. 
Consider the interval [A^, + !]• We have two possibilities: 1) R{N + 1) > 
R*{N). In this case, R*{N + 1){N + 1)^ = R*{N){N + 1)3 > R*{N)N^. 2) 
R{N + 1) < R*{N). In this case, R*{N + 1)(A^ + 1)^ = R{N + 1)(A^ + 1)^ > 
rIn)N^ > R*{N)Ni. □ 

For r > 0, consider the strip 

U{r) := {z e C:\Re z\ <r} . (4.32) 

Generally, we will restrict our attention to small magnetic fields, that is /i e [/q : = 
U{ho) where ^0 will be taken small enough. For instance, /iq < ]^ so that the 
results on the restricted phases can be used in Uq. 
We define the domain of analyticity for a contour: 

Ur:=UiR*iViT)))nUo, (4.33) 

and for a chain X: 

Ux - n . (4.34) 

That is, Ux = Urmii^, where F™'^'' e X has the largest interior. Notice that the 
domains Ur, Ux depend on 9. Set V{X) := ^(F"^'^^) = max{V(T) : F G X}. The 
main result of this section is the following. 

Proposition 4.1. Let 6 E (0, 1), e > 0, c > small enough. There exists 
Pi = l3i{d,e) such that for all (3 > f3i, the following holds. For each chain X, 
h (-^ uj^{X) is analytic in Ux- Moreover, 

\\u^{X)\\u, < uoiX) , II A^+(X)||^_^ < uoiX) , (4.35) 

where uJo{X) is defined via the weights Ao(-P) and Po{T) given in ()4.37|) - ()4.38|) 
hereafter, and satisfies ()4.21|) . 
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Before starting the proof of Proposition l4.H we give explicitly the weights Xo{P) 
and poiX)- These weights are defined such that they can be used throughout the 
section, also when bounding the first derivative of uj^{X). As will be seen, the 
non-trivial part of ^''"(r) will be bounded by: 

e-(intr; +ar] 



r.^r r. ^ <e'3''"e5iri. (4.36) 
0+(intr;-ar) " ~ 

Using dmH), ||pi(r)||c7o < e-^li^lle^^^ol'^leil^l, this suggests to define the weight 
Po(r) in the following way: 

Po(r) := Di/?|r|^e-(^-^)^lirile2^'^«irieir| . (4.37) 

d 

The term Di/3|r|''-i has been added to take into account other contributions, 
especially when studying the first derivative. For clusters we get, using the defi- 
nition of A+(P) and ^W^, 

||A+(P)||h+ < lk+(P)||H+ell-'-(^)l'-+ 

< |^o(i^)|el""(^)l < IMne'"- < D2\MP)\ = UP) ■ (4.38) 

The numerical constants Di, D2 are assumed to be fixed and sufficiently large, in 
order to cover all the cases that will appear in the sequel. 

Lemma 4.3. Let 9 G (0, 1), c > 0, and e > be small enough. Assume 2hQ < 
|(1 — 9)p (p is the Peierls constant). There exists (3i = (3i{9,e) such that for all 
P ^ Pi, the hypothesis ()4.20|) of Lemma \4 . 1\ are satisfied. 

Proof. Define a new weight for polymers (see (j3.52j) ): 

uo{P) := u;o(P)e(^+^(''+^))l^l . (4.39) 
If (3 is large enough, we can proceed as in (j3.6(jp and get 

P90 P90 P90 

This shows the first inequality of ()4.20j) . For the second, we use the Peierls 
condition ||r|| > p|r| fProposition 12. 2|) . This gives 

Po(r)e(^+^)l[ri'l < Di/? |r|^e-(i-^)^''irie2/3'^°irieirie('=+^)l[r]'l (4.40) 

r:[r]i90 r:[r],30 

<Di/3 Y |r|^e-^(i^^)^^irieirie('^+^)l[r]'l. (4.41) 
r:[r],90 

Since |[r]/| < 3'^|r|, a standard Peierls estimate allows to bound this sum by | as 
soon as (3 is large enough. □ 

Until now we have denoted by = er(/3) the small function appearing in the 
study of the restricted phases. Similarly, we denote by ec = ec(/?) the small 
function appearing in the study of chains. These two parameters are assumed to 
have a common bound maxje,., ed < e, which is small. 
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Consider the weight i^'''(r) given ()4.12|) . We can use the hnking procedure for 
the partition functions 0^(intr; =Fc"r); yielding 



(4.42) 



Proof of Proposition \4.1\ The proof will be done by induction. We say a contour 
r is of class n if V{T) = n. A chain is of class n if V{X) = n. 
Consider a contour F of small class (say, of class smaller than Z"'). Then the last 
ratio appearing in ()4.42|) equals 1. We bound uj~^{T) at h = x + iy E Ur- First, 



-2phv{r)\ 



< e 



2(3\x\V{r) 



< e 



2(iR*{v{r))v{r) 



< e 



2/3i?(y(r))y(r) 



< e 



(4.43) 



where we used the definition of the isoperimetric constants K[-) given in (|2.47|) . 
Then, write 



^r(Pintr(+^r))fe ^ ^r(^intr(+^r))h 2r{Vintri+(^r)) iy 2r{'PLri-(^r))iy 
MKtA-^r))h M'P,^tA+^^))^y MKtri-^r)U MKtri-^T))h 



(4.44) 



The middle term has modulus 1 by symmetry (see ()3.38|l ). The two other terms 
can be treated as follows: 



log 



^r(Pi;tr(+^r)) 



ly 



s+iy 



We used Theorem 13.11 Proceeding as in ()4.43p . we get 



ZriKtA-crr)) 



liri 



Ur 



< |x|e^y(F) , 

(4.45) 

(4.46) 



when P is large enough. Altogether this gives 



Ur 



< 



I Pi 



(F)||f;,e^^lir|'eiiri < e-(^^emn^2(,ho\r\^2l\T\ ^ ^^(p) ^ ^4^47) 



Since ||A~''(P)||f/Q < Xo{P), we have shown the first inequality of (j4.35p for chains 
of small class. For the derivative, a Cauchy estimate (any disc centered at h E Uq 
with radius ^ is contained in H^) gives 



|-A+(P)||^^,<16||AnP)b,. 



(4.48) 



For contours 
d 



.-(r) i-C^rm) i-,E^ - 2mr) + ^ log 



Using y(F) < |F 



d 
d-1 



(this is a consequence of Lemma I2.10|) and 
^r)) 



d i,g.^^(^rntr( 



< 2e,y(F) , 



(4.49) 
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this gives the upper bound 



Ur 



<6/5|r|^||a;+(r)||c/,, 



which imphes, as can be seen easily, that 



(4.50) 



(4.51) 



With Lemma [4. H this shows the proposition for chains of small class. Suppose it 
has been shown for chains of class < n. By this induction hypothesis, 1)4.211) and 
Lemma [B. 11 a cluster expansion can be used for the partition functions containing 
chains. Let X be a chain of class n + 1, and consider F G X. The treatment of 
the restricted phases is the same, and we must study the ratio 

^K,r{+ar))h S(A'.;,r(+ar)), E{Xr^,^{+crr)).y E{X+^{-ar))i 



ny 



H(A'+r(-^r))/. S(A'i-,r(+ar))., H(A'+r(-ar))., S(A'+r(-^r))h 



(4.52) 



Again the middle term has modulus 1 and the rest is treated using the induction 
hypothesis. 



l^gS(A'r^r(+^r)). 



This implies 



d 

ds— logS(A'i^tr( 







ds 



(yr))s+ry < \x\e,V{T) . (4.53) 



< 63 



For the weight of F, we thus have (compare with ()4.47p ): 

||^+(r)||f/, < e-(i-e)/3||rile2/3ho|r|g3i|r| ^ ^^^^^ 

For the derivative, use again the induction hypothesis, and bound 

:+^r)) 



dh 



(4.54) 



(4.55) 



(4.56) 



It is easy to check that ()4.5U|) still holds which, in turn, implies (|4.51|) . This shows 
the proposition. □ 

4.4. Pure Phases. In the last section we gave for each chain X a domain Ux 
in which the weight uj~^{X) behaves analytically. The size of the domain Ux 
shrinks to {Re /i = 0} when the size of the largest contour of X increases. In the 
present section we show that the weights uj~^{X) can actually be controlled when 
< Re /i < where is fixed, independently of the size of X. This treatment 
is standard and was first introduced by Zahradm'k |Zlj . 

We consider only chains of type +, the case — being similar by symmetry. Define 



U+ := {z e C : < Reh < h+} , (4.57) 

where < /i+ < min{|^, |} is fixed (p is the Peierls constant). In Sectional 
domains will have to be made optimal, with 9 close to 1, but here we choose 
6* := |. The main result of this section is the following 
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Proposition 4.2. Let e, c > be small enough. There exists (32 = /?2(e) such 
that for all (3 > (32, the following holds. For each chain X of type +, h ^ uj^{X) 
is analytic in U^, and 



\\u;^{X)\\u^<uJo{X), 



(4.58) 



where uo{X) satisfies ()4.21|1 . 



Proof. Since f/+ C H^, clusters P and restricted pliases are under control. For 
each r, we use the representation ()4.12j) (rather than ()4.42|l ). The main ingredient 
of the proof is the following lemma, whose proof is standard and can be found, 
e.g. in |Zlj or |FPj (with minor modifications due to the fact that we are working 
with analytic restricted phases rather than ground states). 

Lemma 4.4. Let (3 be large enough. Then for each contour T of type +, we have 
0"''(intr; — (Tr) on and 

e-(intr;+ar 



< 63 



2|r| 



C/4 



e+(intr;-(Tr) 

The proof of Proposition 14.21 finishes by using Lemma f4. II 



(4.59) 
□ 



5. Derivatives of the Pressure 

In this section we prove Theorem 11.31 adapting the mechanism used by S.N. 
Isakov for the Ising model. Although estimates of Theorem 11.31 hold for the pres- 
sure density p^, we will always work in a finite volume A, and obtain bounds on 
the derivatives of the pressure that are uniform in the volume. As in the preced- 
ing section, we assume 7 G (0,70) is fixed. 

We consider a box A = [-M, +MYnZ'^, with M large, chosen so that A G C^'^. 
Outside A we fix the spins to the value -|-1, i.e. we consider the set fi^, defined 
in (|4.2p and the associated partition function Z^{A) defined in (j4.H) . The finite 
volume pressure p^ is defined by 

<,A^=;^iog^^(^)' (5-1) 

Clearly, this function equals the density pressure of p.l6|) in the thermodynamic 
limit. Consider the set C~^{A) of all external contours of type + associated to 
the set Qj^. Remember that V{r) = |intr|, where intP denotes the union of all 
components of with label — . The family C~^{A) can be totally ordered, with 
an order relation denoted ^, such that V^(F') < V^(F) when F' ^ F. When F is 
not the smallest contour we denote its predecessor (w.r.t. ^) by i(F). 
For a given external contour F G C~^{A), consider the set 

fi^(F) := {o"A G : F' ^ F for all external contour F' of crA+A<=} , (5-2) 

and define the partition function 

e+(A):= Yl exp ( - /?ifA(aA+Ac)) . (5.3) 
crAen+(r) 
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When r is the largest contour then clearly Op(A) = Z^{A) and when F is the 
smallest contour, we define 6^pj(A) := Z+(A). We also introduce the following 
set in which the presence of F is forced: 

f2^[F] := {(Ja G : r' ^ r for all external contour F' of (7\+\c (5.4) 
and F is a contour of cta+ac} . (5.5) 

The partition function 6p;](A) is defined as ()5.3p . with ^^aI-'^] place of r2^(F). 
We have the following fundamental identity: 

e+(A) = e+,)(A) + e+(A). (5.6) 

A crucial idea of Isakov is to consider the following identity. 

zH^) = z-iA) n #7XT- (5-7) 

reC+(A) ^m^^^) 
Then, the logarithm is written as a finite sum: 

logZ+(A)=logZ,+ (A)+ ^ u+{T), (5.8) 

r6C+{A) 



where 



<(r):=log|£^. (5.9) 

Using (|5.6p we can write wj^iT) = log(l + (^^(F)), where 

e+ (A) 

Non-analyticity of the pressure is examined by studying high order derivatives of 
the functions fXi'^) at /i = 0, using the Cauchy formula 



To obtain bounds on (y9^(F)'''^''(0), we exponentiate fXO^) and use a stationary 
phase analysis to estimate the integral. The contour C will be chosen in a k- 
dependent way. If the domain t/p ^ in which ^fXC^) analytic is too small, 
then no information (not even the sign!) can be given about ipX{Tf\0). 

For a while, consider the structure of the partition function Gp^j(A). We write 
A = extAF U F U intF, where extAF := extF n A. By construction, extAF and intF 
are at distance at least / > 2R. We will therefore consider extAF and intF as 
independent systems (see Figure Ej). The sums over configurations on extAF and 
intF can be done separately, yielding 

e+ (A) = p(F)e+r)(extAF; +Acar)e-(intF; ar) , (5.12) 
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A 



+ 



r 




p(r) 



I I 



(extAr; or) 



Figure 6: The decomposition H5.12|l of the partition function (A). 



All the contours of these partition functions are at distance larger than / from F, 
and have an interior smaller than V{r). The point is that we control these func- 
tions for h E Ur, where [/r C C is a domain that depends only on the volume of T. 

The program for the rest of the section is the following. In Section IKTTl we show 
that ^\{r) can be exponentiated, using the results of Section HI We then use a 
stationary phase analysis and obtain upper and lower bounds on some derivatives 
of fXi^) "^aO^) at /i = 0. In Section IH^ we fix k and take the box A large 
enough. For a class of contours called /c-large and thin, the A;-th derivative of 
u~l(T) can be estimated from below, using the results of Section I^TTl This gives 
a lower bound on p^^'*(0). In Section 1331 we show that for p^y^' operations 
limA and {■)^'^^''~{0) commute, leading to the proof of our main results. 

5.1. Study of the Functions (y9^(F). The following lemma requires the main 
results of Sections El and El After that, the proof of non-analyticity of the pressure 
will essentially follow the argument of Isakov (see [l^, |I2] or |FPj ). 

Lemma 5.1. Let 9 G (0, 1), (3 large enough. Then the following holds. For all 
contour F e C~^{-^) with V{r) ^ there exists a map h ^ (7^(F)(/i) analytic in 
the strip Ur, such that for all h E Ur, V2^(r) can be exponentiated: 



^+(F) = exp ( - P\\T\\ - 2phV{T) + 2f3V{T)gliT)) . 



(5.13) 



Moreover, we have the following local estimate 




(5.14) 



The functions 5i are such that limp /-oo Si = 0. 




(5.15) 
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Proof. Consider 0|^j(A). We have seen how to re-sum over configurations on 
extAF and intF. We write 

+ e+r)(extAr;+AcCTr)e+(intr;-ar)e-(intr;+fTr) 

^^^^^ = p^^^ ej;^) enintr;-.p)- ^'-''^ 

All the volume contributions coming from the first quotient will be shown to 
vanish. The partition functions 0T'^p^(extA; err) and 6^ (intF; =FC"r) are of the type 
(j4.5p . We can therefore apply the linking procedure and obtain a representation 
of the form ()4.18|) for each of them: 

ejr) (extAF; +^.ar) = e'^''l^^*'^^IZ,(P+ur)S(^etur) (5-17) 

e±(intr; Tar) = e^^''''^^^ ZAV^^,^)E{X,^,^) , (5.18) 

where we omitted, in the notation, to mention that the families of polymers 
and chains always depend on the boundary conditions specified by +ac and ar- 
Moreover, the family Af^^^.^p contains chains X that satisfy V{X) < V{r). In the 
same way: 

e+r)(A) = e^'^I^IZ.(P+)H(A'+), (5.19) 

where the families Vj^ and depend only on the boundary condition +ac. 
Using the definition of p(T), it is easy to see that ifXi^) form ()5.13|) . 

where g\(X) is defined by 

2f3V{r)g+{r) := -pY,^iM^)-m^\ + ^ogQr + logQ,, (5.20) 
where w(o"j) = —hai, and the quotients Qr, Qc are defined by 

^ ru\ ■_ ^^•('^cxtAr)^r(^i^tr) ^r(^intr) (^ot\ 

z^ivt) zAvl^r) ^'-''^ 



-(^e ttAr)-(^intr) S(A'.; 
S(A'/) S('^intr; 



Since all the families of chains involved contain contours with an interior smaller 
than r, /i I— i> gXiJ^) is analytic in the strip Uy (by Proposition 14. 1|) . Rearranging 
the terms of the cluster expansions for Q,. leads to 

ioga=iog|:j^+ ^'■(^)+ E ^^(^)- E ^^(^) 

pn[r]n¥^$ Pn[r]H^0 Pn[r]fl^0 

Notice that the volume contributions from extAF cancel, and that the three sums 
are boundary terms. By symmetry, the quotient equals 1 at /i = 0, and so 

|logg.(0)| <3e,|[r]^|. (5.23) 

For the derivative, using ()3.63|) gives 

d 



, <2e,\/(r) + 3e,|[r]^|. (5.24) 
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The same computations can be done for Qc- Clusters of chains are denoted X. 
The contributions from extAF also cancel. Indeed, consider the difference 

uj^{X) - u^X) (5.25) 

Using Lemma 1231 there exists for all X\ G '^ettA? "with (i(Xi,r) > R, a cluster 
X2 e X^, X2 n extAr ^ 0, d{X2,T) > R, such that io+{Xi) = uj+{X2)- We are 
thus left with 



xn[r]fl^0 xn[r]H^0 xn[r]«^0 

Using symmetry, 

|loggc(0)| < 3e,|[r]«|. (5.26) 
For the derivative, a similar treatment gives 



<2ecV{r) + 3e,\[rU\. (5.27) 



Estimates (jO^ and yield 

2/3y(r)|(7X(r)(0)| < 3(e, + e,)|[r]^| < (5i(/3)/3||r|| (5.28) 

where 5i(/3) := 3'^"''^/3~^(er + ee)p~^ (p is the Peierls constant). We get ()5.15|) by 
setting 62{f3) := /^-^(e, + ee). □ 

We are now in position of computing derivatives of the functions (/?^(r). The 
main ingredient is the following theorem, which is a consequence of Corollarv lA.il 
of Appendix ^ The first proof of this result was given by Isakov; it is nothing 
but a stationary phase analysis applied to the Cauchy integral giving the A;-th 
derivative at /i = of a function of the type e~'^^~^^-^^^\ 

Theorem 5.1. Let r > 0, F{z) = exp(— C2; + hf{z)) where I < b < c, and f is 
analytic in a disc {\z\ < r}, taking real values on the real line, with a uniformly 
bounded derivative: 

snp\f'iz)\<A<^. (5.29) 

There exists k^ = kQ{A) such that the following holds: define /c+ = r(c — 2b\/~A). 
For all integer k G [A;o, A;+] there exists G (0, r) and > satisfying 

k k 3 I 1 , , 

<rk< TT > 77^ rr, < Cfc < —= , (5.30) 



c + bA c — bA 10 j2'Kcrk y/cTk 

such that 

= ^ [ = A:!^F(-r,) . (5.31) 

In particular, F^''\0) > 0. Moreover, if f satisfies the local condition 

bf{0)<-arc, (5.32) 
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with a G (log 2, 1), then for all k G [ko, k^] and A sufficiently small, 

(log(l + F))^'^(0) = il + a-e--^^')F^'\0), (5.33) 
where a is a bounded function of k, c, b and C = C('^) > 0. 

In Lemma we have put (fiXO^) form e"'^^"'"*-^*^^). In order to satisfy ()5.29|) . 

we must introduce a distinction among the contours. Consider the function 62{P) 
of (UTTT^ . 

Definition 5.1. A contour T G C+(A) is thin z/ |r| < ^v^(r), and fat z/ 
noi thin. 

Now, any thin contour F satisfies, when /3 is large enough, 

A.,tir)l^<2W).Ai3)<2.. (5.34) 

Lemma 5.2. T/iere exzsis /cq swc/i that when (3 is sufficiently large, the following 
holds. For all thin contour T, define 

k+{T) := 2(3V{T)R*{V{T)){1 - 2^1) . (5.35) 

Then for all integer /c G [/cq, ^+(r)]; have 

{-l)''uX{rf \0) > ^(2/?\/(r)D_) V(^+^^(^))l'^l' (5.36) 
(-l)'=n+(r)('^)(0) < 20(2/?y(r)D+) V(^-''i(^»ll^ll , (5.37) 
where lim/3_,oo D± = 1. 

Proof. Let F be a thin contour. Consider </'^(F) in its exponentiated form (|5.13p . 
We apply Theorem lO with c = b= 2f]V{T), f = ^|(F) - r = R*{V{T)), 

and A = A{P). ()5.34p guarantees (j5.29p . There exists = rfc(F) and = Cfc(F) 
such that 

(-l)V+(F)('=)(0) = fc!^V^+(F)(-r,). 



(5.38) 



Using the analyticity of ^'^(r) in Ur, we have with ()5.30p 

(^+(r)(-rfc) = e-^lirile'=^^e'=^^(r)(°)e'=(^'^(r)(-'-^)-^i(^)(°» (5.39) 

>e-/^lirileiTAe-'5i/^l!rile-i^'= (5.40) 

= e-(i+5i)/3||r|lg/Cg"T^^_ (5_41) 

Using the Stirling formula and the estimates for r^, Ck, we get 

(-l)Vl(r)^'>) > ^(2/3y(F)Z}_)'e-(i+'^^)'3||r|| ^ ^5^42) 



where 

2A 

D_{p) = {1- A)e~^ . (5.43) 
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Using ()5.14|) we can satisfy ()5.32|) : 

6/(0) = 2(3V{T)gl{Tm - /5||r|| < -(1 - 6^)(3\\T\\ (5-44) 

< -(1 - 6i)2pV{T)R*{V{T)) (5.45) 

= -(l-(5i)rc. (5.46) 

In (|5.45|) we used 

liril > -—^——V(T)^ > 2V(T) > 2V(T)R*(V(T)) (5.47) 

We can thus use fl5.33|) once /? is large enough. This gives the lower bound ()5.36|) . 
The upper bound is obtained similarly. □ 

5.2. Derivatives of the Pressure in a Finite Volume. In this section, we fix 
k large enough. When a thin contour satisfies [A;o,A;+(r)] 3 k then u^(r)('')(0) 
can be estimated with Lemma 15.21 To characterize this class of contours, we 
introduce a fc-dependent notion of size. 

Definition 5.2. Let A; G N, e' > small enough. A contour T is A;-large if 

V(T) > Vo{k) where 

V,ik) := f^^Millfik)^^ , (5.48) 

where K{oo) was defined in Lemma \2. 91 T is k-small ifV{r) < Vo(/c). 
Let No{e') be such that for all > iVo(e') (see LemmaR3|l. 

^ ^ r < R*{N) < r • (5.49) 



(1 + e') 2K{oo)Nd 2K{oo)m 

Let /c_ = k_{e', 7) be such that when k > k_ then Vo(^) ^ ^o(f')- This definition 
implies that when k > k_, we have for all /c-large contour F 

k4T) = 2mm - 2VA)R*{V{T)) > ^.^/^~f^V (r)^ > k . (5.50) 

K (ooj(l + £ ) 

That is, the k-th derivative of a A;-large thin contour can be studied with Lemma 
15.21 The dependence of k_ on 7 comes from the bound K{oo) > c_7. We 
therefore have lim^\o k^ = +00. 

Proposition 5.1. Let 9 be close to 1, (3 large enough. There exists a constant 
Ci > and an unbounded increasing sequence of integers ki,k2, ■ ■ ■ such that for 
large N, we have whenever A is sufficiently large, 

^ ' ^ 4(r) > (CiK(oo)^/5"^)'"/t;v!^ (5.51) 



lAI dh''^ ^ 
' ' rGC+(A) 



/i=0 



Proof. Fix e > small and consider the sequence (F7v)Ar>i of Lemma l2^ We 
have limjv^oo ^(r^v) = +00 and when is large enough. 



d-l 
d 



[1 - e)K{oo) < ' < (1 + e)K{oo) . (5.52) 

I 
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The sequence {kN)i^>i is defined such that the contribution from the contour Fjv 
to P^A^''(0) is close to maximal. Let 



Cat 



d-1 



ir 



N\ 



(5.53) 



Since lim7v-»oo ^(Tat) = +00, we have lim7v-+oo ^Af = +00. From now on we 
consider large enough so that ()5.52|) and ()5.57|) hold and k]^ > maxj/co, A;_}. 
When considering the /cTv-th derivative, we use the following decomposition: 

E = E 

rec+{A) 



(5.54) 



r6C+{A) rec+{A) rGC+{A) 

fcjv— large, thin fcjv— small, thin fat 



We show that the dominant term comes from F^r, which belongs to the first sum, 
and that the two other sums are negligible. To see that F^r appears in the first 
sum, we first show that F^r is A;7v-large. Indeed, if 9 is close to 1 and e, e',y4(/3) 
are small, 



Vo{kN) < 



< 



K{oo){l + e')d-l 
6(1 -2VA) d ^' 
1 l + e' d-l 



l^fA) 1 



Then we show that F^v is thin: 



IF 



N\ 



< 



1 IIF 



A^ 



< 



d 



1 



d 
d-l 



(5.55) 
(5.56) 

Finally, we assume A is large enough in order to contain at least ||A| translates 
of Ftv. Then we apply Lemma to u|^(Fjv). Using ()5.52p . 



Id z 



y(F^)'=^e-(i+^^)^lir^ll > ((1 



> lil-€)K(oo)^^^kN 



d-l 13 



e)K(cx))||Fyv||j"-^ e 

g (l+5i)^(fciv+l) 

d 



id 



kNld-l 



where c(/cAr) > C^kj^^ and we used the Stirling formula. Since 

(-1)^^^M+(F)('=^)(0) > 

for all k^-large thin contour, we can bound the first sum from below using only 
the contributions coming from the translates of Fjy. We get 



(5.58) 
(5.59) 

(5.60) 
(5.61) 



E 



r6C+(A) 

fcjv— large, thin 



> 



h=0 



c{k 



20 



2'=^K(oo)^''^/3"^''^ 



d 



'1 - e)e~^'D^ 



(5.62) 
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Consider now a A;Ar-small thin contour, i.e. R*{V{r)) > R*{Vo{ki\i)). Using the 
Cauchy formula with a disc of radius i?*(Vo(A;Ar)) centered at /i = 0, 

^ ^'"||n+(r)||^,. (5.63) 



~R*{Vo{kN))- 

Lemma 5.3. Setting ai = ai{9, (3) := p^\l - 6(1 + A{(3)) - 6i{(3)). If (3 is large 
enough, we have ai > and the bound 



ll<(r)lkp < 

Proof. Using (pn!?|) . and dHISl, 



e-/5"i|r| 
_ e-/3°i|r| 



|¥'l(r)||t/, < sup e-/^(i-^i)lirile2/3(i+^)|i^-'^l^(r) < e-^^^l^l < 1 



heUr 

where we used the definition of the radius of analyticity: 



9. 



sup \h\V{T) < R*{V{T))V{T) < RiViT))V{T) < -\\T\ 

hdUv ^ 

The proof finishes by using the Taylor expansion of log(l + x). 



(5.64) 
(5.65) 

(5.66) 
□ 



A standard Peierls estimate implies, when (5 is large, the existence of a constant 



Ca such that 



^ e-^"^|r| < C4IAI . 
rec+(A) 



Using the Stirling formula, it easy to see that kN^.k^ 1 ^ ^ 1. , 
contribution from the A;7v-small contours is then bounded by 



(5.67) 



< kMl^e—"''. The 



1 



rec+(A) 

fcjv— small, thin 



< 



h=0 



ed-i 



9 



2VA 



_d_ 

kNld-i (5.6^ 



Since -At > e<*-i, the comparison of the square brackets of ()5.68|1 with those of 



()5.62j) shows that if 9 is close to 1, if e, e' are small, and if P is large enough, then 
the contribution from the /cAr-small contours is negligible in comparison to the 
A;Ar-large ones. 

We are then left with the contribution of the fat contours. We can use a Cauchy 
bound 



:<(r) 



h=0 



< k\ 

< k\ 

< k\ 



R*{V{T)) 
/2K(l)x'= 



V" 



V{T) 



<(r)llc/r 

e-/3°i|r| 



2K(1) 
9 



-Pai\V\ 

e-/5«i|r| 
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Then a Peierls estimate leads to 

J2 |r|^e-"^^iri < |A|5^LSe-"i^^ < \A\{a[f^)--^T{^ + l) , (5.69) 
rec+{A) L>i 

where r(a:) is the Gamma-function. Using the Stirhng formula, it is then easy to 
show that the contribution from the fat contours is bounded by 

^ J2 4(r) _ <{K{l)p--^D{k)Yk\^^, (5.70) 



lAI 



dh>' ^ ' h=0 

reC+(A) 

fat 



where lim^^oo D{k) = 0. The fat contours can thus always be ignored. This 
finishes the proof of Proposition 15.11 □ 

With ()3.64j) . we get the lower bound, for a large enough box A, 

b^A'^O)! > (CiK(oo)^r^)'"A;^!^ - C^-k^l (5.71) 

> (C_7^/5^^)''^A;^!^ - C^^'kNl • (5.72) 

We used the lower bound K{oo) > c_7 from Lemma 12.91 Notice that we could 
extract the contribution of the translates of Fjv to P^a^''(0) without knowing 
its explicit shape. This is where our formulation of the isoperimetric problems 
differs from the one of Isakov. Notice also that the lower bound (I5.72j) shows how 
non-analyticity is detected in finite volumes. 

5.3. Thermodynamic Limit; Proofs of Theorems (TT^ and [1.31 To extend 
the bounds we have on p1^^^^\o) to the infinite volume limit, we first show that 
in the strip [/+ the derivatives of the pressure are uniformly bounded. 

Lemma 5.4. Let (3 he large enough. There exists C+ > such that for all k > 2, 
sup \\p^^A^\\u+ < (C+7^r^)'A:!^ + C^kl . (5.73) 

A 

Proof. Like in Section we can fix ^ := |. The term C^kl comes from ()3.(i4|l . 
Consider "UaI^) and the representation ()5.1fij) of (F). From Lemma f5. If (pXC^) 
is analytic in Ur- From Proposition 14.21 and Lemma f4. 41 it is also analytic in U+, 
i.e. in U+ U Ur- Proceeding like in the proof of Lemma f5.H we get 

e,^P)(extAF; o-r)6"^(intF; -cxr) 



0;r)(A) 



< sup e-/^R«'^irie3(^'-+^^)l[r]«l (5.74) 
u+ heu+ 

^ ^3{tr+ecW]R\ _ (5.75) 

Assume 3'^+^{er + < |. Using 

\H{r)\\u^ < e"/5||rile'3'^+iriel^l < e""^^!^! < 1 . (5.76) 

Notice that unlike in ()5.65p . a2 in independent of 9. This implies that W^lT) is 
also analytic in [/+ U Ur- Set = min{ai, 02}. Using a disc of radius R*{V{T)) 
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around each /i e [/+, we have 

IHirY'^Wu, < V^(r)||f/,uf/, (5.77) 

<«(!m)K(r).^±^ (5.78) 

<,,(?^)V|-^^ (5.79) 

We used the isoperimetric inequahty of Lemma r2.10[ Remember that K{1) < C4.7 
(Lemma 12. 9p . and that / = The proof finishes like for the upper bound on 

fat contours. □ 

Corollary 5.1. For all h' eU+U {Re h = 0} and for all k eN, 

= hm p^f{h') = lim p^J:\h) . (5.80) 

AyZ'^ h\h' ' 

Proof. We show fl5.8()|l for A; = 1. By definition, 

#-(/zO = lim^2i^^^±4^M^ (5.81) 

^ <5\0 6 

= lim lim ^ i ^^'^^ ^ (5.82) 

= lim^hm (p+«(/z') + ^ptA^(M5))'^) , (5.83) 

where lim^^o h{6) = h'. The following lemma will allow to permute the limits 
lim^^o and limy^^^-zd. 

Lemma 5.5. Let, for all N eN, 6 > 0, biy{6) = + ciy{6) , such that \cn{S)\ < 
DS uniformly in N, and limAr_^oo ^Af('^) = b{S) exists. Then lim m-, 00 a n and 
lim^^o K^) exist and are equal. 

Proof. We first show that lim^^^o K^) exists. Let (Sk) be any sequence (5^ > 
such that limfc^oo 5fc = 0. Then we have 

\b{5k) - b{6k')\ = I lim (c^(4) - c^(4'))l < Di^k + 4') (5.84) 

and so lim^^oo b{6k) exists. Fix e > 0. There exists N^^g such that if > A^^^^^ 
then \biy{S) — b{6)\ < e. We then have 

b{6) - e - D5 < liminf qat < lim sup cat < b{5) +t + D6, (5.85) 

which finishes the proof, once we take e 0, 5 — > 0. □ 

Using the fact that the second derivative is uniformly bounded on (Lemma 
15. 4|) . this shows the first equahty in ()5.80p . For the second, we only need to 
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consider the case where h' = 0. 



pW-(0) = hm^^2^^1^^ (5.86) 



5 5 



hm 

(limip«(M5)))+ip«-(0), (5. 



(5.87) 



where h{5) G [|, 5] and hm^^o h{6) = 0. This shows 



P 



(0) = \im £\h{5)), (5.89) 



(1), 

7 v-y ^^^^r-^ 



which extends easily to any sequence h \ 0, since derivatives of any order are 
uniformly bounded on ?7+. □ 

We can then complete the proofs of our main results. 

Proof of Theorem \l.!^ - The bounds on p^^^^i^) of ()5.72|) and Lemma EiH extend to 
the thermodynamic limit using Corollary 15.11 □ 

Proof of Theorem \l.^ - Using the symmetry p^{h) = p^{—h), we can write, for 
m > 0, 

f-y{m) = sup [hm — p-yih)) . (5.90) 

h>0 

By the Theorem of Lee and Yang, h i— s> ^^(/i) is analytic in {Reh > 0}. If m* : = 
p^^'^{0), then : (0, +oo) (m*, 1) is one-to-one, and for all m G (m*, 1), 

f-^/{m) = h{m)m — p^{h{m)) , (5.91) 

where h{m) is the unique solution of the equation piy\h) = m. The GHS in- 
equality (see |GHSj ) states that for all /i > 0, 

(0) >pf (/i) > 0. (5.92) 



P 



(2). 
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Since p^\h) ^ for all h > 0; the biholomorphic mapping theorem implies 
that m I— *• h{m) is analytic in a complex neighbourhood of each m G [m* , 1). So 
/y, which is a composition of analytic maps, is analytic on (m*, 1). 
Since p^^'^{0) = \imh\o p^\h) exists (Corollary 15. 1|) . we can extend p^^ to 

[0, -|-oo) by setting p^\o) := p^^''~{0) and we can extend h{-) to [m*, 1) by 
setting h{m*) := 0. 

We now show that has no analytic continuation at m* . Assume this is wrong. 



"'^"'^Let (7 : D ^ C be an analytic function, zq G be a point such that ^'(zo) 7^ 0. Then there 
exists a domain V <Z D containing zq- such that the following holds: V' = giV) is a domain, 
and the map g : V V' has an inverse : V ^ V which is analytic, and which satisfies, 
for all uj € V', g^^'(uj) = (g' {g~^ (uj))) . The proof of this result can be found in |Rem| . pp. 
281-282. 
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and denote by h{-) the analytic continuation oi h{-) at m*. Using the definition 
of h{m), we compute 

7n\, , him) — him*) 

h^^Um*) = hm ^ — ^ 5.93 

m\m* m — m* 

h(m) - him*) 

= hm ^ ^ -^T^ 5.94 



So /i(-) maps a neighbourhood of m* on a neighbourhood of 0, and hS^\m*) ^ 0. 
The biholomorphic mapping theorem imphes that h{-) can be inverted in an open 
neighbourhood of /i = and that this inverse, which coincides with p>l^ on the 
positive real line, is analytic, a contradiction with the non-analyticity of at 
/i = (Theorem mni). □ 

6. Conclusion 

Our analysis has lead to the following representation of the pressure for /i > 0: 

P,{h)=pt,^{h) + s1:^{h), (6.1) 

where is the restricted pressure. As we have seen in Section pj!"^, which 
describes a homogeneous phase with positive magnetization, behaves analytically 
at /i = 0. On the other side, contains the contributions from droplets (con- 
tours) of any possible sizes, and is responsible for the non-analytic behaviour of 
the pressure at /i = 0. Non-analyticity can be detected only in the very high or- 
der derivatives of s+, although sl^ contributes essentially nothing to the pressure 
when 7 is small. Indeed, can be expressed as a sum over clusters of chains, 
and each chain contains at least one contour. Since the length |r| of a contour is 
bounded below by the size of a cube C^'^ , we have 

\\s^,\W<ae-'^^<-\ (6.2) 

where a, 6 > are constants. 

For the pressure, the Lebowitz-Penrose Theorem takes the form (see |Prj ) : 

Po{h) ■=\imp^{h) = inf ^^(m) , (6.3) 

7\0 mG[-l,+l] 



where 



(phifn) '■= —hm — —I{m) . (6.4) 

2 p 



The bound ()6.2|) implies, for /i > 0, 

Poih) = limp+ (/i) = inf (j)h{m) . (6.5) 

From this last expression, the analytic continuation of the pressure, in the van 
der Waals limit, at h = 0, can be understood easily: for /i > 0, 0/i(m) has a 
unique global minima at m*{h,P) > 0. When h < this minima is only local. 
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h = Q 



Figure 7: The function (j)h{m) for different values of h. For h > 0, the vertical line 
represents the value taken by po{h). When /i < it represents its analytic continuation, 
given by \im^\^op+^. 



but provides the analytic continuation at /i = 0. The identity shows that 
the constraint on the local magnetization, in jo^^, has the effect of selecting the 
minima m*(/i,/3), which is global when h > and local when h < 0. This 
mechanism is illustrated on Figure [7| When 7 > 0, this scenario breaks down: 
droplets are well defined, and they are all stable at /i = 0, creating arbitrarily 
large fractions of the — phase. As we saw, this gives a contribution k\'^ to the 
A;-th derivative of the pressure. 

Appendix A. A Stationary Phase Analysis 

The following theorem is a generalization of a result due to Isakov jll]. Let 
Vp{t) := {z e C : \z - t\ < p}. 

Theorem A.l. Let p > 0, F{z) = exp{—cz + bf{z)) where 1 < b < c, and f is 
analytic in a disc Vp{0), with a uniformly bounded derivative: 

sup \f{z)\<A<^. (A.l) 

There exists ko = ko{A) such that the following holds: let t G 'Cp(O) and define 
k^ = {p — \t\){c — 2b\/A). For all integer k E [/cq, k^] there exists z^ = r^e^'^'' G 
T>p{0) and c^ G C such that 

F^'\t) = k\j^Fi-z, + t). (A.2) 

When Imt = and f{z) takes real values for real z, then ^pk = and Imck = 0, 
and we have the estimates 

3 1 1 1 

< Recfe < ——= and llmcfcl < . (A. 3) 



10 V27rcr^ y/cr^ .Jcr\ 

I \ bA k cos k cos ipk 

tan(/?fc < — and — - < ru< — - . (A. 4) 

' ^ I - c - c^bA - c-bA ^ ' 

We have not indicated, for notational convenience, the dependence of rk,ifk,Ck 
on t. A consequence of this Theorem, for t G (— p, +p), is given after the proof 
in CoroUarv I A. 11 Our theorem improves significantly the original result of |TT| . 



55 



since we show that derivatives of the function can be obtained anywhere in the 
disc of analyticity Vp{0), although we don't use it in our proof. In the course of 
the proof, we make clear the fact that the stationary point = rfcC**^* is solution 
of a system of equations (see ()A.9|) - ()A.10|) ). whereas Isakov considered only the 
point t = 0, and there a single equation suffices to find Zk since ipk = 0. Since 
this result is at the core of the proof of non- analyticity, we have explicited every 
step of the proof. 

Proof of Theorem \A.l\ We use the Cauchy formula. Define k G (0, 1) by /^p = 
For all r G (0, p(l — k)) we have 

k\ 2txi JsvAt) - ^) 27rr'= Jq e^'fi 

where we have used the parametrization z := t + re*'^ for &Dr{t). Our aim is 
to extract the main contribution to this last integral. The integrand, because of 
the form of F, has a maximal value for close to tt. We thus make a change of 
variable, (f':=ip — n,to obtain 

where 

0(r, (/?) := cr cos Lp + 6Re/(— re**^ + t) 

'(/'(r, ip) := cr sin p + blmf{—re^'^ + t) — kp . 

If t G M and / is real at real points, f{z) = f{z) and therefore lmf{z + t) = 
—lmf{z + t). By symmetry we get F^''\t) G M. The core of the proof is to choose 
r in a specific manner. This is a standard stationary phase analysis. To this end, 
we will need an estimate on the second derivative of /. Using Cauchy allows to 
obtain, for all Zq G T>p{0): 

/ (^o) = — / 7 r^du = — / — d^ , (A.6) 

where r' > is such that I^qI + < P- Using the uniform bound \f'\ < A we 
obtain (we optimize taking the largest possible r', namely p — \zo\) 

ir(^o)|<— (A.7) 
P ~ I ^0 1 

Now, set t G T>p{0), \t\ = up, and consider the map p t-^ /(— re"^ + 1). A direct 
computation yields 

-ire''^f'{-re''^ + t) 

re'^f{-re'^ + t) - (re*^) V"(-re*^ + t) . 



d 

dp 
d^ 



dp 



f{-re'^ + t) 
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Using (jA.7|) gives the bound 

d2 



sup 



dip- 



:f{-re'^ + t) 



<rA + r' 



A 



rA 



p{l — k) — r 1 



(A.8) 



We now turn to the existence of a saddle point. 



Lemma A.l. Let t G Vp{0), \t\ = up. Then for all k G [0, p(l - k)(c - 2b\/~A)], 
the system 



^0(r,^) = O 
d 



dip 



ip{r, V?) = 



has a solution (r^, ip^) with and pi^ satisfying the following estimates: 

I , bA k cos pk k cos pk 

tan Pk < — CLiT'd — —TT' "^fk"^ 



c + bA 



c-bA 



c-bA 

Proof We exphcit (ICTl . (IXlOll : 

sin (/?(c — bRe/'(— 2 + t)) — cos blmf'{—z + t) = 0; 
r cosp[c — bRef'{—z + t)) + r sin p blmf'{—z + t) = k . 

These two equations are equivalent to 

ksinp = rblmf'{—z + t)) ; 

k cos p = r{c — bRef'{—z + t)) . 



(A.9) 
(A.IO) 

(A.ll) 



(A.12) 
(A.13) 



Then, we see that any solution of the system ()A.9|) . ()A.10|1 . satisfies ()A.11|) . To 
show that there exists a solution, we first solve ()A.13|) locally for some fixed p> G 



lip _ 



(— |, +|) (so that cos (p > 0). Define the map r i— ^(r, p) := r(c — 6Re/'(— re 
t)). Since / is analytic, its real and imaginary parts are C°° with respect to r > 
and p (see jRemj), so ^ is C°°. We have ^(0, p) = 0, and 

e(p(l - k), v?) = p(l - «:) (c - bRef'i - p{l - K)e^^ + t)) 

> p{l- K){c-bA) 

> p{l - k){c - 2by/A) >k>kcosp 

which proves the existence of some G (0, p(l — k)] such that ^(r^, p) = k cos p. 
Notice that we also have that 



k cos p 



c-2b\fA 
< —— < 1 . 



p(1-k) p{l - k){c - bRef) - c - bA 



(A. 14) 



12n 



The definition of fc+, witli 2\/ A instead of y A, ensures the strict inequality < 1 
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We can then show that the solution is unique, by verifying that ■^^{r,ip) is 
strictly positive at r^. First, 

— e(r, f) = c- bRefi-re'^P + t)- r6— Re/'(-re^'^ + t) 
or or 



see 



(US) 



c - bRef'{-re'^ + t) + r6Re(e^^/"(-re^^ + t)) 
bA 



> c 



1 - 



At r = r^p we get (see fjA.14p ) 
d 



dr 



cy/A 8c 

> — > 



2-VA- 9 



which proves uniqueness of r^. The continuity of i— > is a consequence of 
the implicit function theorem. We turn to the second equation, and set r = r^. 
Using again equations ()A.12|) . ()A.13|) . we have 



tan(/? 



On (— |, |) the function 



blmf'{r^e'^) 
c - bRef {r^&'fi) ' 

blmf'{r^e'^) 
c-bRef'{r^e'^) ' 



-bA bA 



is continuous and takes its values in the interval ( 
exists a solution {rk,(pk), fk '■= r^^, of ()A.12|) and ()A.13|) . 



c-bA' c~bA 



). Therefore there 



Notice that we have explicit bounds on ipk, such as 



simpkl < \ tanipk\ < — , jcosv^fc] > 



4 



TT 



□ 



(A.15) 



and that we can estimate, at r = (see flA.Sjl and ()A.14|l ) 

d2 



sup 



dip- 



< 



fkA , C — bA 5 c r-r 

<rkA-:^^ — _< rfcVI. (A.16) 



p{l-K) 



Ib^fA - M ~ 9 6 



We now examine ()A.5|) when r = rk- Defining Zk = rkC^'^'^, we extract the value 
taken by the integrand at Zk'. 



F{-Zk + t) 1 



2n 



{''fe,¥')-<}!>(r-fe,<^>fe)+i(i/'(rfc,(/3)-V'(rfc,i/3fe)) 



We will estimate the integral in Ck by decomposing [— tt, +7r] in two parts. The 
first is [-7r,+^]\[-f,+f]. 



13 



We have not yet shown that this solution is unique. 
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Lemma A. 2. For all 6 > 0, there exists ki = ki{6) such that for all k > ki we 

have 



1 

2^ 



+ 



dip 



< 6- 



1 



V2ncrk 



(A.18) 



Proof. We have \e'^'('(n,v)-Hrk,Vk))\ = i. pirst, consider the interval [f,7r]. On 
this interval, cosy? < y{ip) where ip t— > yi^p) ■= cos ^ — sin|((/? — ^) (we have 
y(7r) = —0, 95 ■ ■ ■ > —1). We can thus compute 

crfc (cos ip — cos ipk) = cr k{cos ip — cos — + cos — — cos ipk) 



- — ^^^'^'-'^ crfc (cos - 

V2 , 2.3 

< crAu) cru , 

2 a' 25 



cos tpk) 



where we used ()A.15|) in the last step. For the other part containing /, 



6(Re/(- 



-re 



lip 



t) - Refi-re'"^" + t)) < b{ip - tpk) sup 

< brkA{(p - ipk) 



—Ref{-r,e'^'^ 
dip 



TT, 2 

< crkA{(p - ^) + ^crfc 



The first part of the integral can thus be bounded by: 



t) 



e 25 



2n 



V2 

e ■ 2 



crfc(¥3-f)+crfcA(<^-f ) 



d(p < 



e 25 



27r 



e '(^"^^""dx 



< 



e 25 



v/2^(^ 



2^) V2?cffc 



< 



2 ySvrcf^ 

once A; is large enough, since cr^ > y+aI^ (^^^ (lA.llj) ). The same can be done 
on [— TT, —j], on which we use the function y{ip) := cos ^ + sin |((y9 — |). □ 



On the interval [— f , +f ], we use Taylor expansions for and -0, around 99 = (/j^. 
We have (r = is fixed) 

1 d^ 

(l){ip) = </)((^fc) + + -{(p - ipkf 

1 d^ 

^{ip) = ^l){ipk) + + —{ip - ipk)^ —0 
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where 93 and </? are both functions of ip. On the interval [— f , +f], we have the 
estimates 

10 5 
CTk < -r^0 < — crfc (A. 19) 



Indeed, since 



d^ d^ 

' -crfc cos + 6— Re /(-rfcC^'^ + f), (A.20) 



we use ()A.16|1 and find 

dV 5 10 



> -cn - ^^crk > -^CTk , (A.21^ 



d(/?2'- - 9 9 

and the upper bound 

dV 7r 5 5 



, ^^<— crfccos — I — y AcTh < — cr^ . (A. 22) 

We can thus compute some upper bound on the integral over [— f , + f ] in as 
follows: 

1 1 r + 00 T f- 

271 J . ^ - 27r ^ ^ 2 9 ' ^ 



00 
1 



The upper bounds on Rec^ and Imc^ can be obtained by taking, say, 5 = | in 
Lemma [a .21 which gives 

\ck\ < <5-L + ^ ^ • 

The lower bound on Rec^ is obtained by dividing [— f,+f] = /i U /2, where 
i\ = [fk ~ li '^k + 1], and 7 = 7fc G [■"!!+!] is defined by two conditions: first, 
fix 7 small enough such that 

2A 

sup I sin Lp\ < . (A.24) 

fell 1 — A 

The existence of such a 7 is guaranteed by (jA.lljl . This first choice implies that 

1 d^ 
sup -tp{(fk)\ < 7T7^sup (crk\smip\ + b| -— Im /(-rfcC*'^ + t) 1) 

<pe/i 2 ^e/i dip! 

( see (HH) < ^7'(crfcY^ + ^crkVA) 

1 2 AT/2VA 5. 

= -TcrkvAi + -) 

2' H-A 9^ 

< -crfc a/A7^ < Y5crfc7^ 
Then, the second condition on 7 is the following: 

^crfc72 = ^ . (A.25) 
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Here, we might have to take k large enough to make sure that 7 G [— f)+f] 
Then, we have a lower bound: 

271" J I, ~ 27r Jj^ 



cos — f 1 10 

see > ^ exp ( - - ■ ycrfc((^ - (pfc)')d<^ 



:(2<|.(J— ^)-l)- 



VSO^ 'V 27 ' ^V2^F^ 
47 1 

> — ^ (A.26) 

~ 100 v/to^ ^ ^ 

The upper bound on /2 = [— ^, ipk — 7] U [ipk + 7, + f ] is obtained easily: 

/ + / )e^(-)-^(-^)d(^ < 2 ■ / ( - ^ ^cr,^2)dx 

27r J_. 27rJ^ 2 9 



< 2 ■ — / exp • -cTkX )ax 

- 2n ^^29^ 

6 . , /so",. 1 

= — (1 - — ^)) , 

< — — (A.27) 



Taking 5 := ^ in Lemma lA. 21 and using ()A.26|) . ()A.27|) gives 

Recfc> )^^ = -^^, (A.28) 

''-UOO 100^^2^ lOv/2^' ^ ' 

which completes the proof. □ 

Corollary A.l. Let p > 0, F{z) = exp(— + bf{z)) where 1 < b < c, and f 
is analytic in a disc Vp{0), taking real values on the real line, with a uniformly 
bounded derivative: 

sup \f{z)\<A<^. (A.29) 

There exists k^ = kQ{A) such that the following holds: let t G (— p, +p) and define 
k+ = {p — \t\){c — 2b\l~A). For all integer k E [ko,k+] there exists > and 
Ck > such that 



F('\t)=k\-^F{-r, + t). (A.30) 



We have the estimates 



3 1 Ik k , 

<c.<^, —i-r<rk<——-. (A.31) 



10 y/2'Kcrk y/ci\ c + bA ~ ~ c-bA 

In particular, (— l)'^F*^'^^(t) > 0. Moreover, if f satisfies the local condition 

-ct + bf{t) < -ape, (A.32) 
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with a G (log2, 1), then there exists function a = a{k,c,b), sup |a| < oo, and 
7 = 7(a) > such that for all k G [k^, k^] and A sufficiently small: 



( log(l + F)f\t) = (1 + a ■ e-"^')F^''\t) 



(A.33) 



Proof. The first part of tlie proof follows from Theorem lA.il For the second part, 
notice that ()A.32|) implies < e~"'"^ < 1. The continuity of / implies that 

in some neighborhood V C T>p{0), V 3 t, we have sup^gy 1-^(^)1 < 
In V we can thus use the Taylor series for log(l + F): 



^apc 



< 1. 



log {1 + F)=J2 ^^^F" 

n>l 



n 



n>2 



■1) 



n+1 



n 



Since the series converges absolutely and uniformly in V , we can derive it term- 
wise with respect to k. We then need to show that the following holds: 

Lemma A. 3. Let a G (log 2, 1). There exists a positive constant Kq < 00 such 
that for all A G (^^, 1) and for all n > 2, k G [ko, k+], 



(A.34) 



where 7 is given by 7 := aX — log 2. The constant A in ()A.29j) has to be taken 
small enough (depending on the value of a). 

Suppose for a while that the Lemma has been shown; we have the following 
estimate: 



n>2 



-1 



n 



<Y^\iF-Y'^\t)\ 



n>2 



< Koe-2''\F^''\t)\J2' 



-Q(l-A)(n-l)fc 



n>2 



< Kn 



-a(l-A)feo 



□ 



' \ — g-a(l-A)fco 

which proves ()A.33|) . 
Proof o f Lemma \A . SI The point is that F is of exponential type. We have 

= g-'^"^ + ^"/(^) = ^-CnZ + b„f{z) ^ (A. 35) 



with c„ = on, bn = bn. For each n = 1, 2, . . . , we can apply Theorem lA.il there 
exists for all /c G [/cq, ^n.+J, where kn,+ = nk+, some rn,k and Cn,k such that 



Cn,k 



-F{-rn,k + tY 



(A.36) 



Notice that [/cq, 3 [kQ,kn^+\ for all n. The constant r„^fc is a solution of the 
equation k = r{cn — bn^ie f'{—r + t)) and satisfies 



k 

Cn ~l~ bnA 



< rn,k < 



k 



br,A 



(A.37) 
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The constant Cn^k satisfies 

3 



< Cnh < 



1 



10 y/2jw~r~k ^/Cnrn,k 
We can then consider, for all k G [k^, k^]: 



(A.38) 



(A.39) 



Notice that when n increases, r„_fc \ and oo. Using ()A.37j) and ()A.38|) . 

we find 



ri,k . 1 + A 
— — < n- 



'^n.k 



I- A 



and 



We must estimate 



cifc 3 V 1 - ^ 



exp {c{nrn,k-ri,k) - ct{n - 1) 



(A.40) 



(A.41) 



+ K^/(-^n,fc + t)- f{-n,k + t))) (A.42) 



Using the definition of rn,k gives 
nrn,k - n^k = k 
= k- 



_c - bRe f'{-rn,k + t) c - 6Re /'(-ri,^ + t) _ 

b(Re f (-r„,fc + t) - Re f (-ri,^ + t)) 
c - 6Re /'(-r„,fe + t)){c - 6Re f'{-n,k + t)) 



2hA 



< 



k 2A 



{c-bAy - c{i-Ay 

We then compute the term involving / (we use twice — f{x')\ < A\x — x'\): 

nf{-rn,k + t)-f{-n,k + t) (A.43) 
= (n - + n(/(-r„,, + t) - /(t)) - (/(-n,, + t) - f{t)) 

< {n-l)f{t)+nArn,k + An,k 

<(n-l)/(t) + - 



c 1 - A 



We thus have 



F{-rn,k + ^ ^,(A)fcg(-rf+fe/(i))(n-l) 



(A.44) 



i^(-ri,fc + t) 

where e{A) = 2A{2 — A){1 — Since pc > k^ > k, we can use assumption 

|02|l and get, for all A G (0, 1), 

^(-ct+bf(t))(n~l) ^ g-afc(n-l) ^ g-aAfc(n-l)g-a(l-A)fc(n-l) (A.45) 
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Since logn — log 2 < ^{n — 2) for all n > 1 we can compute the following bound 
sup n'^e-"'^^^"-^) = sup e'=(i°g"-"^('^-i)) (A.46) 

n>2 n>2 

< supe'=('°*5 2-'+°^+"(^-"^)) (A.47) 

n>2 

< gfc(log2-aA) ^ g-Cfc^ (A.48) 

where we used the fact that A is chosen such that C = > 0. Putting our 
bounds together we bound ()A.39|) . when A is small, by 

f l±^e^(A)^ g-|Cfcg-|Cfcg-a(l-A)fc{n-l) ^ ^^g-iCfcg-a(l-A)fc(n-l) _ 



A 



□ 



Appendix B. Cluster Expansion 



Consider a countable set V whose elements are called animals, and denoted 
7 G P. To each animal 7 is associated a finite subset of Z'^, called the support of 
7. Usually we also denote the support by 7. In the cases we consider, the support 
is always an /^-connected set. Assume we are given a symmetric binary relation 
on V, denoted ~. We say two animals 7, 7' are compatible if 7 ~ 7'. When 7 and 
7' are not compatible we write 7 7^ 7'. We assume that the following condition 
is sufficient to characterize incompatibility: for each each animal 7, there exists 
a set 6(7) C Z'^ such that if 7 7^ 7', then ^(7) n b{y) ^ 0. 

To each animal 7 G P we associate a complex weight G C. The partition 
function is defined by 

^(^)^= E n (B-i) 

{-i}CV 76(7} 
compat. 

where the sum extends over all sub-families of T> of pairwise compatible animals 
(we assume this sum exists, which is the case in every concrete situation). When 
{7} = 0, we define the product over 7 as equal to 1. We are interested in studying 
the logarithm of the partition function. To this end, we define the family D of 
all maps 7 : P ^ {0, 1, 2, . . . }. The support of 7 is the set {7 G P : 7(7) > 1}. 
Usually we also denote the support of 7 by 7. We will also write 7 9 x if the 
support of 7 contains an animal whose support contains x. A map 7 G P is a 
cluster of animals if its support can't be decomposed into a disjoint union 5*1 U 5*2 
such that each 71 G S\ is compatible with each 72 G 6*2. Formally, the logarithm 
of the partition function has the form (see e.g [FT] ) 

logS(P) = J]a;(7), (B.2) 

where the weight of 7 equals 

c.(7)=a^(7)n^(7r^'^ (B.3) 

7G© 
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The functions a'^(7) are purely combinatorial factors. They equal zero if 7 is 
not a cluster. The following is the technical lemma that gives explicit conditions 
for the convergence of the development ()B.2|1 . The proof is standard and can be 
adapted from |PTj . 

Lemma B.l. Let uJo{'j) be a positive weight such that 



where < e < 1. Define ujq{'j) as in ()B.3|1 with u!q{'j) in place of uj{^). Then 
there exists a function ri{e), \im^_Qr]{e) = such that 



Typically, in the cases we consider, the weights are maps z 1— u^j^z), an- 
alytic in a domain ^ C C, and there exists a positive weight <-iJo(7) such that 
||'^(7;-)IU — 1^0(7) for all 7. Lemma FB.!! thus implies that the series ()B.2jl is 
normally convergent on A. This guarantees analyticity of the logarithm of S(X'), 
by a standard Theorem of Weierstrass (see e.g. |Remj ) . 
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